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FORBUORD 


The space shuttle main engine, a reusable space propulsion syst^, is the 
most sophisticated rocket engine in use today. However, on the premise that 
the United States will depend on the shuttle and its derivatives as its prin- 
cipal Barth- to-orbit transportation system for the next several decades, the 
National Aeronautics and Space Administration (NASA) has decided to continually 
upgrade SSKB capabilities. Therefore, NASA has initiated a research and tech- 
nology program for advanced high-pressure oxygen-hydrogen rocket technology. 
This program is intend‘d to establish the basic discipline technology and to 
develop new analytical tools for an orderly evolution of reusable rocket 
syst^BS . 

A substantial porcion of this technology program is devotea to the struc- 
tural integrity and durability aspects of reusable space propulsion systems. 
This effort %rill i^rove the understanding of aerothermodynamic loads, struc- 
tural dynamics, and fatigue, fracture, and life issues in the SSME and its 
derivative engines. The overall approach is to examine and modify existing 
analytical methods, to develop new methods where needed, and to verify them in 
systoaatic test programs. The research is conducted through industrial con- 
tracts, university grants, and in-house activities at the Lewis Research 
Center. 

To provide a forum for representatives from government, industry, and 
acadoaia to learn about and discuss the latest findings and progress toward 
ii^roved structural integrity and durability of reusable space propulsion sys- 
tems, a two-day conference was held at Lewis in June 1985. This publication 
contains extended abstracts of the papers presented at that conference. 

Stanley J. llarsik 
Conference Chairman 
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OAST TBCHNOL06Y PROGRAM - DURABILITY EFFORT 


F.W. Stephenson, Jr. 
RASA Headquarters 
Washington , 0 . C . 20546 


The NASA Office of Aeronautics and Space Technology has in place a program 
directed toward technology advances needed to improve service life, perfor- 
mance, and operations of reusable rocket engines such as the space shuttle main 
engine (SSHE), its derivatives, and advanced engines of this class. The pro- 
gram is focused on the development of (1) accurate analytical models for 
describing flow fields, aerothermodynamic loads, structural response rotor- 
dynamics, bearing characteristics, etc., from which life predictions codes can 
be evolved; (2) the technology for longer life components such as bearings, 
seals, turbine blades, and combustors; and (3) advanced instrumentation to 
better measure the engine internal operating environments and to develop diag- 
nostic sensors for real-time condition monitoring of critical engine compo- 
nents. The program is supported by a strong effort in materials and processes 
development and evaluation. 

An integral part of the reusable rocket engine technology program is the 
structural integrity and durability effort, which is the subject of this con- 
ference and which is focused on the hot-gas flow system of staged combustion 
engines like the SSMB. This presentation provides a brief overview of the 
overall technology program and the relationship between the overall program and 
the structural integrity and durability effort. 
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OVERVIEW OF AEROTHBRMODYNAMIC LOADS DEFIEITIOfir STUDY 
Louis A. Povinelli 

Rational Aeronautics and Space Administration 
Lewis Research Center 
Cleveland 9 Ohio 44135 


The objective of the Aero thermodynamic Loads Definition Study* is to 
develop methods to more accurately predict the operating environment in the 
space shuttle main engine (SSME) components. Development of steady and time- 
dependent, three-dimensional viscous computer codes as well as experimental 
verification and engine diagnostic testing are considered to be essential in 
achieving that objective. The steady, nonsteady, and transient operating loads 
must all be well defined in order to accurately predict powerhead life. 
Improverents in the structural durability of the SSME turbine drive systems 
will depend on our knowledge of the aerothermodynamic behavior of the flow 
through the prebumer, turbine, turnaround duct, gas manifold, and injector 
post regions. 

The study was begun in October 1983. The various elements in the study 
were described in December 1983 at the initial meeting of the SSME Durability 
Activity at Lewis. During the first year, emphasis was placed on the computa- 
tions of the turbine losses, the redistribution of the radial temperature pro- 
file through the turbines, and the flow behavior in the turnaround duct. A 
viscous computation was also carried out for the fuel turbine rotor, and work 
was begun on a time-dependent turbomachinery code. Experimentally, an unsteady 
heat transfer apparatus was designed, and instrumentation was fabricated for 
installation in a low-aspect-ratio turbine. Highlights of the first year’s 
accomplishments were as follows: 

(1) Computations of the fuel turbine losses indicated a potential improve- 
ment of 8 to 10 points by reducing blade roughness and platform leakage. 

(2) Oxidizer turbine loss could be reduced 5 to 7 points by redesigning 
for correct rotor incidence angles. 

(3) Turbine inlet temperature profiles were not significantly changed 
through the turbines. It was concluded that the high solidity led to rela- 
tively lightly loaded surfaces and reduced secondary flow generation. Hence 
temperature redistribution was reduced. 

(4) Viscous computations of the first-stage rotor profile, using flow 
conditions supplied by Rocketdyne, revealed a flow separation on the pressure 
side. Subsequent computations, using newly computed stator exit conditions, 
did not rc^veal any separation. 

(5) The first manned orbital flight (FMOF) and full power level (FPL) 
turnaround ducts were analyzed by using a viscous computation. After a 
detailed study of the pressure losses, flow separation, diffusion, and 


*This study is part of the activity of the Fluid and Gas Dynamic Working 
Group of the SSME Technology Program. 
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acceleration through the turn, a new geometry was computationally defined 
within the available space envelope. The modified design of the turnaround 
duct had minimum pressure loss and reduced the potential for flow separation. 

(6) Substantial effort was devoted to fabricating and installing high* 
response heat flux gauges on a low-aspect~ratio turbine. This apparatus will 
be used to determine the response of the blade surfaces to temperature varia- 
tions introduced by upstream disturbances. 

The results obtained during the first year of the study were presented in 
two workshops (March and November 1984) and at the Advanced High Pressure 
O 2 /H 2 Technology Conference (June 1984), all held at the Marshall Space 
Flight Center. A summary of the results was also presented at the SAE Aero- 
space Congress Meeting (November 1984). A list of references is included. 

The material to be presented in this session will update these results. It 
represents the effort since October 1984 and is focused on the turbine compo- 
nents. It is anticipated that the upcoming yearns effort will be directed 
toward the unsteady aspects of the turbine flows. 
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SSME FUEL PREBURNSR TWO-DINENSIOHAL ANALYSIS 


Thomas J. VanOverbake 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland. Ohio 44135 


The durability of the SSME turbine is strongly affected by the temperacure 
profile leaving the prebumer. Experimental measurements of this ^^emperature 
profile are difficult to obtain because of the severe operating environment 
within the SSME. Therefore a designer attempting to tailor the flow entering 
the prebumer to conform to acceptable temperatures must either engage in an 
expensive and difficult experimental effor or be overly conservative in the 
design. With recent advances in computatic .al capabilities, a more promising 
design approach is to use a reacting flow computer model to predict the turbine 
inlet temperature profile. This is the objective of the current study. A 
series of calculations wero made, using a reacting flow code, to assess the 
sensitivity of the turbine inlet temperature profile to variations in the flow 
entering the SSME prebumer. 

A two-dimensional teaching type of combustion code was modified to work 
with prebumer geometry. The program solves the fully elliptic, steady-state 
equations of motion. Either hybrid or bounded skew upwind differencing can be 
used. The pressure distribution was found by the pressure- imp lie it split 
operator (PISO) predictor-corrector technique. Closure for the equations of 
motion was provided by using the turbulent kinetic energy k and energy dis- 
sipation to obtain an effective eddy viscosity. The density was calculated 
from the ideal-gas law. Combustion was modeled by using the eddy breakup model 
of Magnussen and Hjertager. This model relates fuel combustion to turbulence. 
Combustion depends on the mixing of species and the dissipation of energy from 
hot combustion products. Combustion of the fuel in the code occurs as a nega- 
tive source term f^L the fuel fraction depending on the dissipation of turbu- 
lent eddies (mixing) and the concentration of fuel, oxygen, ana product in each 
grid cell. 

The prebumer was modeled as a scries of coaxial nozzles injecting into a 
chamber. 'he walls were approximated by a fine rectangular mesh. Thus the 
rear wall of the prechamber model is a series of steps over which the gases 
flow to the powerhead turbine. A 70-by-53 grid was used. The density of the 
radial grid points was increased for the igniter and for the exit flow to 
ensure an accurate solution. The turbine inlet temperature profile was taken 
at 0.33 m, but calculations were performed to 0.5 m to ensure that exit bound- 
ary conditions for the code would not affect the temperature profile. The 
cooling flow along the outside wall, which eventually enters the prechamber, 
was included in the nozzle flow nearest the outside wall. 

For most of the calculations the preburner inlet was divided into five 
regions of premixed reactants. The igniter was modeled as an inlet with a 
mixture of hydrogen and oxygen at low temperatures or a mixture of hydrogen and 
water at 1000 K. Changing the inlet condition prodvxed little change in .he 
turbine inlet temperature profile oven though the flow was greatly changed in 
the inlet computational area. Also, changing the equivalence ratio of the 
inlet corresponding to the igniter did not change the turbine inlet temperature 
profile to any appreciable degree. One reason is that the igniter takes up the 
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smallest area of the inlet flow and cont.ibutes the smallest part of the flow. 
Changing larger areas of the flow than that of the igniter did change the tem- 
perature profile. Decreasing the equivalence ratio of any of the other inlet 
conditions produced temperature increases in the flow corresponding to that 
portion of the total inlet. Changing the inlet condition for the largest area 
had the greatest effect on the turbine inlet temperature profile. Also» 
increasing equivalence ratios decreased temperatures: the fuel-side SSME pre- 

burner operates so fuel rich that increasing the fuel flow decreases tempera- 
tures and decreasing the fuel flow increases tei«q)eratures . Decreasing just the 
outer inlet equivalence ratio resulted in a series of inlet boundary conditions 
with decreasing equivalence ratio toward the outer prebumer wall. This run 
generated a temperature profile with temperatures increasing from the inner 
wall to the outer wall. A comparison calculation was made with increasing 
equivalence ratios toward the outer wall, keeping the overall mf.ss flow simi- 
lar, The resulting temperature profile was skewed in the opposite direction 
from the previous calculation. 

The inlet turbulence was increased to 30 percent from 5.5 percent and 
decreased to 0.55 percent to see if the inlet turbulence had an effect on the 
temperature profile. Calculated temperatures differed only for a few grid 
points downstream of the nozzle inlets, and the turbine inlet temperature pro- 
file was essentially the same. The code was also run with five burned inlet 
conditions. That is, the corabustion had already occurred, the reaction was not 
needed, and the temperature was a function of mixing. The turbine inlet tem- 
perature profile was essentially unchanged, even though inlet flow conditions 
were vastly different. 

Finally, the number of inlet boundary conditions was increased to repre- 
sent a coaxial series of alternating pure hydrogen and oxygen nozzles along 
with the igniter. The turbine inlet temperature profile was changed but only 
for a few of the outer grid points. Most of the turbine inlet temperature 
profile was unchanged even though temperatures varied widely downstream of the 
inlet nozzles because of the extreme variation in local equivalence ratios. 
Since the actual inlet for the prebumer is a collection of more than 200 noz- 
zles rather than a small number of coaxial nozzles, the calculation with 
alternating pure aydrogen or oxygen coaxial nozzles corresponds to a worst 
case. Thus the actual temperature profile should be closer to that generated 
with the five premixed inlet conditions. 

The temperature profiles generated by these calculations with the eddy 
breakup combustion model show that the turbine inlet temperature profile is 
affected by changing the equivalence ratio (if the affected area is large 
enough) but that other aspects of modeling, such as modeling separate oxygen 
and hydrogen nozzles or modeling the inlet conditions as burned or unbumed, 
have minor effects. 
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COMBUSTION MODEL 


THE PRESENT CODE USES THE EDDY BREAKUP COMBUSTION 
MODEL WHICH RELATES COMBUSTION TO THE DISSIPATION OF 
TURBULENT EDDIES AND THE CONCENTRATION OF SPECIES. 

Hg ♦ 5 Og HgO ♦ 1.21x10® J/kg 

h - mp Hr ♦ Cp^ J 

^p.m* 2 “iSjI 

S.I * fi, "=ll * "=2l^ * ♦ C^lT^ ♦ Cs,T^) 


THE COMBUSTION IS MODELED AS A SOURCE TERM 
FOR THE MASS FRACTION OF FUEL, nip 

Sfflp » -/^ln[ Amp. A 5^^ • AB sxqIC+1^ 

A and B ars reaction conatante 

e/K may be thought oF ae an eddy IIFetIme 

THE RATE OF COMBUSTION IS 



FUEL PREBURf«R IN^^CTOR CONFIGURATION 



THE PREBURNER INJECTOR WAS DIVIDED INTO FIVE ZONES 


ZONE 

NASS FLOW C kg/a) 

VELOCITY 

TENP 

EQUIVALENCE 


°2 

«2 

Cffl/a) 

CK) 

RATIO 

IGNITER 

.209 


77.6 

1000.0 

9.2 

A-B 

4.70 

4.81 

23.9 

162.8 

8.18 

C-D 

8.46 

7.92 

24.7 

162.5 

7.49 

EFG 

19.75 

17.44 

23.3 

162.4 

7.07 

H 

7.52 

8.27 

25.6 

162.9 

8.79 


T„ =165 K T„ -125 K 

”2 °2 

Zona H Inciudaa liner cooling 


POOR QUALJry 
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RADIAL DISTANCE (m) RADIAL DISTANCE (m) 


OE POOR y 


VELOCITY VECTOR DIAGRAM FOR BASEUNE CALCULATION 
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AXIAL DISTANCE (m) 

VELOCITY VECTOR DIAGRAM USING 5 PREMIXED INLET B.C. 




Ofii 



0612 r 


9996 


•0000 ' 




^ ^ rm ^ 


: EE = ====“==“::s^^^^^ 
• » ^ ^ ^ ^ ^ ^ ^ ^ ^ / ' 


; siasjisaaaasasas; 


X 


X 


X 


X 
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AXIAL DISTANCE (m) 

VELOCITY VECTOR DIAGRAM USING 4 UNBURNED INLET B.C. 
THE INLET B.C. FOR ZONE A-B WAS ALSO USED FOR THE IGNITER 


ZONE 
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(MET*:R) 03 , , . ^ (HEUQ'' 


TURBINE INLET TEMPERATURE PROFILE AS A FUNCTION OF INLET B.C. 



OP!S, EQUIVALENCE RATIO TEMPERATURE 

\SBJNE CALCULATION USING 5 PREMIXED INLET B.C. CHANGING THE INLET 
TURBULENCE DIDN'T EFFECT THE TEMPERATURE PROFILE NOR DID 
MODELING OF THE INLET FLOW AS HOT BURNED GASES. 



ORIG. EQUIVALENCE RATIO TEMPERATURE 

CrlANQNG THE INLET CCNDR ION FOR THE IGNITER TO THAT OF ZONE A-B 
DIDhTT CH\NGE THE TEMPERATURE PROFILE NOTICEABLY, WHILE CHANGING 
B.C. FOR THE LARGER INLET ZONES DID CHANGE THE TEMPERATURE PROFILE 
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INLET TURBINE INLET TEMPERA1URE PROFLE AS A FUNCTION OF INLET B.C. 
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ORIG. EOUI^AiLfNCE RATIO TEMPERATURE 

CHANGING THE INLET B.C. FOR ZONE E-F-G TO THAT OF ZONE H EFFECTED THE 
TE'MPERATURE PROFLE THE MOST. ZONE E-F-G IS THE INLET i OR ALMOST 50% 
OF THE TOTAL FLOW. 



CRIG. EQUIVALENCE RATIO TEMPERATURE 

TEMPERATIFIE PROFILES RESULTING FROM A SET OF INCREASING OR 
DECREASInIG EQUIVALENCE RATIOS. THE INLET B.C. WERE CHANGED 
SO THAT THE OVERALL MASS FLOW WAS SIMILA.H 
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TEMPERATURE TEMPERATURE 

BASELINE 5 PREMIXEO INLET B.C. 

ALTERNATING UNMIXED INLET NOZZLES 

THE TURBINE INLET TEMPERATURE PROFILE RESULTInKS FROM A LARGE NUMBER 
OF ALTERNATING PURE HYDROGEN OR OXYGEN INLET NOZZLES WAS ONLY FFFECT 
AT THE TOP. EVEN THOUGH INTERMEDIATE TEMPERATURES DID VARY CONSDERABL 



TEMPERATURE TEMPERATURE 

BASELINE INLET B.C. IGNITER B.C. SANE AS ZONE A-6 


USING LOWER REACTION CONSTANTS IN THE EDDY BREAKUP CONBUSTION 
MODEL RESULTED IN VARIOUS TEMPERATURE PROFILES 
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SSME TURBOPUMP TURBINE COMPUTATIONS* 


Philip C.E. Jorgenson 
Sverdrup Technology , Inc . 
Middleburg Heights, Ohio 44130 


A two-dimensional viscous code developed at the NASA Lewis Research Center 
is being used to predict the flow in the SSME high-pressure turbopump blade 
passages. The rotor viscous code (RVC) employs a four-step Runge-Kutta scheme 
to solve the two-dimensional, tl;in-layer Navier-Stokes equations. The Baldwin- 
Lomax eddy-viscosity model is used for these turbulent flo^ calculations. A 
viable method has been developed such that the relative exit conditions from 
an upstream blade row are used as the inlet conditions to the next blade row. 
The blade loading diagrams are compared with the meridional values obtained 
from an in-house quasi-three-dimensional inviscid code. Periodic boundary 
conditions are imposed on a body-fitted C-grid computed by using the GRAPE 
(GRids about Airfoils using Poisson’s Equation) code. Total pressure, total 
temperature, and flow angle are specified at the inlet. The upstream- running 
Riemann invariant is extrapolated from the interior. Static pressure is spec- 
ified at the fxit such that mass flow is conserved from blade row to blade row, 
and the conservative variables are extrapolated from the interior. For viscous 
flows the no-slip condition is imposed at the wall. The normal momentum equa- 
tion gives the pressure at the wall. The density at the wall is obtained from 
the wail total temperature. 

To compute the flow for a given blade row, inlet conditions were based on 
the mass-averaged relative exit conditions from an upstream blade row. Inlet 
conditions to the first stator were taken from conditions given by Rocketdyne. 
Mass flow was then matched by varying the exit static pressure. This process 
was repeated for each of the blade rows. The calculations were done at midspan 
of each blade. 


DISCUSSION 

The results from the RVC were compared with those previously obtained by 
rising an inviscid quasi-three-dimensional code at the mean span section of the 
blade. *^he comparisons were made for the engine balance flew conditions given 
by Rocketdyne for the SSME high-pressure fuel turbopump. Figures 1 and 2 show 
the computational grids used in compu^'ing the first stage. Figure 3 and 4 show 
the loading diagrams for th.. stator and rotor. Originally a separated region 
was computed on the pressure suricce of the first-stage rotor (figs. 5 and 6). 
In the present calculation the flow haj an adverse pressure gradient but does 
not quite separate (fig. 7). The original cC’^outation did not account for the 
effects of an upstream blade row; therefore the lirst-stage rotor had a lower 
inlet flow angle. The computational grids for the second stage are shown in 
figures 8 and 9. The RVC results of the second- stage stator compare quite well 
with the inviscid code results in figure 10. The difference in the loading 
diagram of the second-stage rotor was due to the different exit static pressure 
needed to match the mass flow of the turbine (fig. 11). Another calculation 
was made to match the exit static pressure of the inviscid code; however, this 
increased the mass flow by 6 percent (fig. 12). Figure 13 compares the 


*Work performed under Lewis task number 82-24-01. 
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velocity diagrams of RVC and the inviscid quasi-three-dimensional code (at the 
mean section) for the SSMB high-pressure fuel turbopump. 

Originally only the two rotors were computed by using the inlet conditions 
given by the inviscid code. The current quasi-coupled results are believed to 
be more realistic. 

Presently, work is being done on the SSME high-pressure oxidizer turbo- 
I»ump. The computational grid and loading diagrams for the first-stage stator 
computed by using the inviscid version of the RVC are presented in figures 14 
and 15. 

A two-dimensional Navier-Stokes code was used to predict the flow in the 
SSME high-pressure fuel turbopump. The relative upstream exit conditions were 
used as input to the successive blade rows, thus coupling the solutions. The 
results obtained by using the Runge-Kutta version of the RVC compared well with 
the results from the inviscid quasi-three-dimensional code. The RVC is good 
for understanding the physics of the flow. Work is continuing on the SSME 
high-pressure oxidizer turbopump. 
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SIMULATION OF MULTISTAGE TURBINE FLOWS* 

M.L. Celestina and R.A. Mulac 
Sverdrup Technology, Inc. 

Middleburg Heights, Ohio 44130 

John J. Adamczyk 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 4413S 


The long-term objective of the numerical simulation of turbine flows is 
to enhance our understanding of the flow phenomena within multistage turbo- 
machinery components. The direct benefit of this activity is improved modeling 
capability, which in turn can be used to improve component efficiency and 
durability. 

The flL f field in multistage turbomachinery components (single stage being 
a limiting case) is extremely complex. The flow field is highly unsteady, with 
time scales ranging from a fraction of shaft speed to several times the highest 
blade passing frequency. The length scales of the flow are also diverse. They 
range from the circumference of the machine to a small fraction of a blade 
chord. 

To assess the difficulty in analyzing the flow field within multistage 
turbomachinery component:., a hierarchy of equations was formulated to serve as 
a base for analyzing these flows. The equation set that provides the most 
complete description is the Navier-Stokes equations. The simplest description 
is given by a set of equations that govern the quasi-one-dimensional flow. As 
one proceeds from the complete description, the number of unkno%eis to be solved 
for increases monotonically above the number of equatiCiis. The development of 
the additional set of equations needed to mathematically close the system of 
equations forms the closure problem associated with that level of description. 
For the Navier-Stokes equation there is no closure problem. For the quasi-one- 
dimensional equation set random flow fluctuations, unsteady (deterministic) 
fluctuations, nonaxisymmetric flow variations, and hub-to-shroud variations on 
the quasi-one-dimensional flow must be accounted for. 

Given the capabilities of today*s supercomputers (including N.A.S.) it 
seems reasonable to attempt to analyze the flow within multistage machinery by 
using a set of equations that describe the time-averaged flow field within a 
typical blade passage of a blade row. With respect to a given blade row the 
flow in a typical passage is steady and spatially periodic from passage to 
passage. There is one such description per blade row. The closure problem 
associated with this description requires accounting for the effects of random, 
unsteady (deterministic), and spatial fluctuations. The development of compu- 
ter codes to solve the typical passage flow equations would enhance our ability 
to simulate the time-averaged, three-dimensional flow field within multistage 
components. These simulations should enhance our understanding of the flow 
physics associated with these devices. For example, a better insight should 


*Work performed under Lewis task number 83-24-01. 
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be obtained into the mechanism by which low-momentum fluid is redistributed 
within a blade row and the effect of this redistribution on blade row perform- 
ance and durability. These solutions should also yield a description of the 
unsteady, distorted flow field encountered by a blade row. 

To date, a hierarchy of equations and their associated closure problems 
has been formulated for use in simulating flows in multistage turbomachinery. 

A computer code has been developed that is based on a Runge-Kutta integration 
scheme suggested by Jameson to solve the inviscid form of the equation govern- 
ing the time-averaged flow in a typical blade passage. The inviscid closure 
problem associated with this equation set is being examined for counterrotating 
propellers, fans, and turbine stages. The preliminary results obtained for a 
counterrotating propeller configuration seem to suggest that rational inviscid 
closure models can be developed for this class of machines. Further work is 
planned in this area on t!.e geometry supplied by General Electric for their 
unducted fan configuration. 

Work will begin soon on including the viscous terms in the code. The 
closure problem associated with this equation set will be examined. It is also 
planned to develop a time-accurate code for solving the Reynolds-averaged 
Navier-Stokes equations by using the Runge-Kutta integration scheme. This code 
will be used to develop data bases for devising closure models for the average 
passage equation set. 
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DIRECT SIMULATION OF UNSTEADY TURBOMACHINERY FLOWS 


SPEED REQUIREMENT BASED ON 15- MIN RUN 
WITH 1985 ALGORITHMS REYNOLDS-AVERAGED 
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THREE-DIMENSIONAL CELL-CENTERED HNITE VOLUME 
EULER SOLVER 


EULER EQUATIONS IN CYLINDRICAL COORDINATES 
(d\Q/dt)dV + S(Xf d4r + XGdAO + XHdAz) - S 
FOR ROTATING SYSTEM 

(aXQ/at)|a6s - (a\Q/dt)\rel - (l{d\Q/dB)\rel 

EULER EQUATIONS FOR POTA '1NG SYSTEMS 
(a\«/a''dV + J^iXFdAr + X(C - mQ)±iO + \HdAz) - S 

JAMESON'S RUNGE-KUTTA INTEGRATION PROCEDURE 

four stage scheme 
second to fourth order cccurocy 
ollowsCFL>1 (Cn. = 2i^) 
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UNSTR4DY FLOW IB MULTISTAGE TURBIBES 


Surya P. Suraiq>udi 
Cleveland State University 
Cleveland • Ohio A4115 


In recent years there has been increasing interest in the development of 
distortion-tolerant engines. The solution to the problem of inlet distortion 
and forced vibration may require major engine redesign and retesting. The 
objective of the present research project is to develop an efficient model for 
the response of a multistage turbine to either a total pressure or total tem- 
perature distortion. 

Each turbine blade row is modeled as an actuator disk. The actuator disk, 
indiich is defined as an artificial device producing sudden discontinuities in 
flow properties, is often used to describe the flow through turbine and com- 
pressor blade rows. The flow approaching the blade row is assumed to be sub- 
sonic and inviscSd. The distortion occurs at upstream infinity and is in the 
form of total pressure or total temperature. Recently M. Celestina, R. Mulac, 
and J. Adamczyk (ref. 1) obtained a solution to unsteady Euler equations by 
using a finite volxime method. With prescribed inlet and boundary conditions 
flow variables such as density, pressure, and velocities can be calculated at 
cell centers in the entire flow region. The inlet and exit of the blade row 
coincide with lines of constant ^ (fig. 1). 

The annular duct code developed in reference 1 is used to obtain the flow 
variables at various cell centers for the region 0 < ^ < n and n ^ 3 < ^ < m 
The characteristic equations for the region n < ^ < n + 3 are obtained from 
equations of continuity momentum and energy. These equations ((1) to (5)) are 
expressed in terms of flow variables such as density, velocity in three direc- 
tions (u, V, w), and pressure p. Nondimensional time and space variables in 
axial (r), radial (r), and circumferential (<p) directions are given by t, 

T)f C- The characteristic equations are expressed in finite difference form 
by equations (7). A multistage Runge-Kutta method is used to integrate these 
equations. Equation (8) gives the details of the time-stepping scheme used in 
the analysis. After carrying out these integrations the relations between the 
flow variables (eqs. (9) to (13)) are obtained. 

By making use of actuator disk theory for the flow across the blade row 
the following assumptions are used in the present analysis: 

(1) Entropy loses are neglected. 

(2) Hub-to-tip radius ratio is high. 

(3) Total enthalpy is constant. 

The exit angle is specified. From these assumptions the relations between flow 
variables at the inlet and exit of the blade row are derived. These relations 
are given by equations ( 1 ^) to (18). Equations (9) to (18) form a set of 
nonlinear equations containing 1C unknown variables (p^^, u^^, p^, Wj^, 

P2> ^2* ^2* P2)* These equations can be further reduced to two 

nonlinear equations ((19) and (20)) containing two unknown variables, u^, U2* 
These two equations can be solved by using a Newton-Raphson method (ref. 3). 

The details of the Kewton-Raphson method are given by equations (21) to (22). 
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A con^uter code based on this analysis is under development to obtain 
flow variables at the inlet and exit of the blade row. The preliminary results 
will be presented at the SSME durability conference. 
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UH8TBADY HBAT TSAVSPBR DDB TO TIMB-DBPimBNT FRBB STREAM VELOCITY 


Rana S.R. Gorla and All Anari 
Claveland Stata univaraity 
Clavaland. Ohio 44115 


There exists a need for a detailed study of the effects of the 
upstream unsteady flow on the heat transfer from a turbomachinery blade. 

One of the most significant studies of laminar boundary layers 
under the Influence of a purely time-dependent, free stream oscillation 
Is done by LIghthlll Cl]. His analysis employed a linearization for 
small oscillation amplitudes. Lin L2D considered the effect of 
finite amplitude oscillation on a flow field. Mori and Tokuda C3D 
Investigated the heat transfer from an oscillating cylinder. 

Recently, Gorla C4D examined the unsteady fluid dynamic characteristics 
of an axlsynwnetrlc stagnation point flow on a circular cylinder 
performing a harmonic motion In its own plane. He presented solutions 
for small and high values of the reduced frequency of oscillation. 

The present work has been undertaken In order to study the 
unsteady combined convection from a horizontal circular cylinder 
to a transverse flow. A coordinate perturbation method Is used to 
transform the governing set of partial differential equations Into a 
system of ordinary differential equations. The free stream time- 
dependent velocity was assumed to be sinusoidal and the boundary 
layer response due to both low as well as high frequencies of oscillation 
will be studied. Currently numerical solutions are being obtained for 
the distribution of the unsteady Nusselt number and the friction factor. 

METHOD OF SOLUTION 


Assuming an Incompressible flow with constant properties and 
negligible dissipation, the governing equations within the framework 
of the boundary layer approximation may be written as 
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In the above equations, x and y represent the distances along the 
streamwise direction and normal direction; u and v the velocity 
components In x and y directions; U© (x,t) the external velocity; t the 
time; v the kinematic viscosity of the fluid; a the diffusivlty and T 
the temperature. 
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The boundary conditions are given by 
y=0: u=v=0 and T = T„ 
y u = Ue (x,t) = Uq(x)*(I + e e^^^) 

T T«, (^) 

The solution strategy Is to subject equations (l)-(3) to the 
linear perturbation analysis, expand the resultant zero and first-order 
equations In a series of x* wh ose coefficients are function of n alone. 

Here, x*=(x/R) and n=y«^2U„/vR. The resultant ordinary differential equations 
are then to be solved numerically. Appropriate unsteady solutions will 
be developed for low frequency and high frequency of oscillations of 
the free stream velocity. 

This work Is currently In progress and It Is hoped to report these 
results In the very near future. 

REFERENCES 

1. LIghthill, M. J. ”The Response of Laminar Skin Friction and Heat 
Transfer to Fluctuations In the Stream Velocity," Proc. Royal 
Society, Series A, Vol. 224, 1954, pp. 1-23. 

2. Lin, C. C. "Motion In the Boundary Layer with a Rapidly Oscillating 
External Flow," Proc. 9th Int. Congress of Applied Mechanics, 

Vol 4, 1957, pp. 155-167. 

3. Mori, Y. and Tokuda, S. "The Effect of Oscillation on Instantaneous 
Local Heat Transfer In Forced Convection from a Cylinder," Proc. 

Third Int. Heat Transfer Conference, Vol. 3, 1966, pp. 49-56. 

4. Gorla, R. S. R. "Unsteady Viscous Flow In the Vicinity of an 
Axisymmetric Stagnation Point on a Circular Cylinder," International 
Journal of Engineering Science, Vol. 17, 1979, pp. 87-93. 


40 



■ N85-27948 

HBH FiU:iLITY TO STUDY UHSTBAOY HAKE EFFECTS ON TURBINE AIRFOIL HEAT TRANSFER 


Robert J. Siuoneau 

National Aeronautics and Space Administration 
Lewis Itesearch Center 
Cleveland, Ohio 44135 


A significant portion of the SSME aerothermal loads program Is directed 
at the heat transfer effects of the unsteady flows, particularly wakes, that 
occur naturally In turbomachinery. Although these phenomena occur In all tur- 
bomachlnes, we feel they will be more severe In the SSME turbines because of 
the high heat transfer associated with the very high Reynolds numbers over the 
SSME airfoils. 


There has been a significant Increase in both analytic and experimental 
research In turbomachinery unsteady flow phenomena In recent years. Coiq>uta- 
tlonal and experimental measurement advances have made this research possible. 
Some very exciting results are beginning to emerge. The next 5 to 10 years 
should see a significant step upward. Most of the experimental work, however, 
will he conducted at Reynolds nttmbers typical of aircraft gas turbines, an 
order of magnitude below the SSME. At the NASA Lewis Research Center we are 
building a research facility to study fundamental unsteady wake effects on heat 
transfer. The facility will operate from SSME Reynolds numbers and span the 
gap down to aircraft turbine Reynol(*s numbers. 


The experimental concept Is shown In figure 1. The test section will be 
a single airfoil In an annular sector with walls contoured to match the air- 
foil. The airfoil will have a chord almost 1 foot long. The wake-forming 
disturbance generator will be a spoked Wheel of small circular cylinders rota- 
ting on a mean diameter of 6 feet. A conceptual layout of the entire rig is 
shown In figure 2. The conceptual design phase Is complete and the final 
design In under way. Construction should begin In the fall of 1985, and the 
rig should be ready for research In the summer of 1986. 


The following Is a brief overview of the expected simulation. Figure 3 
Is a schematic of wake behavior In a turbine stage. When a rotor blade passes 
through the stator wake, the suction surface sees the wake first. Since the 
wake absolute velocity Is less than the mean, the velocity triangle shows that 
the wake relative velocity vector will be deflected around from the design 
stagnation point toward the suction surface (l.e . negative Incidence). A 
careful examination shows exactly the same behavior for the rotor wake on the 
second stator. Thus, Whether the wake Is stationary or moving In space, the 
effect Is the same. Figure 4 represents the simulation. As with the real case 
the bar wake In the simulation Impinges on the suction surface first. The bar 
drag will also deflect the velocity vector toward negative Incidence. The main 
similarity variables in this experiment are Reynolds number ard Strouhal num- 
ber, defined as follows: 


&iVem . 

Re =» ; 


2vfL 
® “ (Vel) 
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Many characteristic velocities and dinensions could be used. Ue selected 
exit relative velocity and axial chord for the passage simulation and inlet 
relative velocity and diameter for the leading edge. The following table shows 
numbers for the SSMB at full power engine balance conditions, for a VASA ver- 
sion of the Energy Efficient Engine (E^), and for the bar wake simulation 
rig. 



SSMB 

e3 

Test rig 

Passage conditions: 

Reynolds number 

7.19x10* 

0.53x10* 

0.77x10* to 7.69x10* 

Strouhal number 

3.81 

2.70 

3.81 

Leading-edge conditions: 

Reynolds number 

0.37x10* 

0.051x10* 

0.062x10* to 0.62x10* 

Strouhal number 

0.59 

1.13 

1.83 

Wakes per passage 

1.6 

1.3 

1.8 


The Strouhal numbers change very little from turbine to turbine While the 
Reynolds nund>ers vary by an order of magnitude. With the test rig the two 
parameters can be varied independently, and the full Reynolds number range can 
be spanned. Finally, in real machines there is always more than one but less 
than two wakes in the passage of one time. This too is simulated. 

In a parallel program being conducted under basic research and technology, 
we have already begun to gather evidence that this approach will yield valuable 
results. Figure 5 shows a rig that is conceptually similar to, but smaller 
than, the new facility. It is a full annular passage approximately 1 foot in 
menn diameter. The test section, rather than an airfoil, is a heated cylinder 
with a splitter plate. Figure 6 shows typical instantaneous and ensenible- 
averaged wakes generated by such a spoked wheel. They look very similar to 
ones reported in the ' Lterature for turbomachinery. Early work in this facil- 
ity indicat<*s that on the average the wakes can be treated as an rms turbulence 
intensity to yield a heat transfer correlation (fig. 7). Although this may be 
useful for design, it does not tell the real story. To see the real effect of 
the wakes, a ceramic cylinder with thin-film sensors painted on the surface 
(fig. 8) was heated in a small radiation furnace and then suddenly injected 
into the flow stream behind the spoked rotor. The resulting instantaneous heat 
flux along the stagnation line of the cylinder (measured directly by the use 
of an analog circuit) Is shown in figure 9. The data indicate a 25-percent 
spike in heat flux at the bar passing frequency. This early work suggests that 
the wake effects on heat transfer are real and measurable and could yield sig- 
nificant high-frequency temperature fluctuations at the surface of a turbine 
rotor blade. The value of this approach to studying the wake-generated 
unsteady flow effects has also been shown in some recent work at the University 
of Oxford in England. 

In summary we feel the rig under design provides an excellent simulation 
of the i&q>ortant unsteady physics in the SSMB turbines. The test section is 
large and accessible, Which will allow detailed measurements to verify the 
complex analyses needed to describe the phenomena. 
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EXPERIMENTAL MEASUREMENTS AND ANALYSIS OF HEAT TRANSFER 

AND GAS DYNAMICS IN A ROTATING TURBINE STAGE 

M.G. Dunn* W.J. Rae, and U.K. George 
CALSPAN-UB Research Center 
Buffalo • New York 14225 


The life-cycle fatigue of a particular component of a high-pressure tur- 
bine stage is influenced by the local heat flux and the resulting metal tem- 
perature excursions experienced by this component during normal operation. 

Flow through a turbine stage is seldom uniform and is often influenced by non- 
uniform incoming gas streams that can result in locally high heat fluxes. In 
the case of the rotor blades passing through the nozzle guide vane wakes and 
passages, these local heat rates translate into local temperature fluctuations 
that will likely affect the fatigue life of the blade. In order to better 
understand the influence of the heat load on the component life, one would 
like to have both the local steady-state heat flux and the high-frequency 
comi-onent of the heat flux. 

Relevant fundamental data of this type can be obtained by making use of a 
full-stage, rotating turbine model currently in use at CALSPAN. The turbine 
being used is the Garrett Turbine Engine Company TFE-731-2 HP stage, and both 
the nozzle guide vanes and the rotor blades are heavily instrumented with thin- 
film heat-flux gauges, as described later. The experimental technique is the 
short-duration, shock-tunnel approach, in which fast-response, thin-lxlm ther- 
mometers are used to measure the surface temperature historias at prescribed 
locations on the turbine component parts. Heat-flux values are then inferred 
from the temperature histories by using standard data reduction procedures. 
Local variations in the stagnation temperature distribution incident on the 
rotor are created by injecting cold gas through the trailing-edge cooling slots 
of the nozzle guide vane. The injected gas (air) is at a known total tempera- 
ture and a known weight flow rate. 

INSTRUMENTATION 

The major emphasis in this measurement program is on acquiring, analyzing, 
and predicting the rotor blade steady-state and time-rasolved heat flux in the 
presence of variations in the incident total temperature distribution. Early 
in the program it was mutually agreed that the blades would be instrumented at 
the 10, 50, and 90 percent spanwise locations with the maximum gauge coverage 
at the 50 percent location. Two blade-contoured leading-edge inserts contain- 
ing four gauges each were fabricated and installed in the midspan region. Five 
additional gauges were installed in the blade tips, and three additional gauges 
were installed on a fourth blade at the 50 percent spanwise location to aid in 
the analysis of the wake-cutting data. The pressure surface instrumentation, 
the suction surface instrumentation, the leading-edge inserts, and the final 
rotor assembly are shown in figures 1 to 4. 

The purpose of this effort was to develop the data interpretation tech- 
niques necessary to interpret the high-frequency heat-flux data. The data 
acquisition system has been configured to acquire phase-locked, high-frequency 
heat-flux data. A previously developed (at CALSPAN) calculational technique 
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w%s integrated with the phase-locked system to enable the calculation of 
instantaneous heat-flux values during wake cutting. These techniques were then 
verified by piggybacking on the Air Force-sponsored Low Aspect Ratio Turbine 
(LART) measurement program. 

The objective of the analytic phase of this research is to provide a 
theoretical/numerical complement to the experimental work. The data being 
gathered under this program are very basic and require continuing examination 
in the light of state-of-art analysis. Thus a key component in assembling the 
CUBRC program in turbomachinery has been the acquisition and exercising of the 
standard NASA flow-field analysis codes (refs. 1 to 3). 

These codes were provided by the NASA Lewis Research Center^ and were 
installed on the CDC Cyber 815 of the SUNY-Buffalo Computer Center. It was 
soon discovered that the storage requirements of the codes strained the limits 
available on this machine, and arrangements were then made to shift to a VAX 
11/780 which had become available for research purposes. Since this is a 
virtual memory machine, no storage complications were encountered. However, a 
small amount of reprogramming was required in order to compensate for the way 
that ENTRY statements are handled by the VAX. These revisions were completed, 
<nd successful test runs were made by the end of the contract period. 

Thus the program objective was completed: the basic NASA flow-field codes 

are now available for use in support of the experimental efforts at CALSPAN. 

CONCLUSIONS 

The measurement program described herein will be run in mid-May 1985. It 
is anticipated that at the time of the June 4-5 meeting, the measurements will 
have been performed. It is doubtful that any significant analysis of the data 
will have been conq>leted at this time, but preliminary comments concerning the 
program will be presented. 
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CALCULATION OF VISCOUS COMPRESSIBLE FLOWS IN SHALL GAS TURBINES 


O.L. Anderson and G.D Power 
United Technologies Research Center 
East Hartford, CT 06108 


The NASA Levis Research Center is conducting a program to develop 
improved technology for the Space Shuttle Hain Engine (SSH£)« This program 
requires not only a component teat program, but also the development of 
accurate and reliable analyses to support the program. Thus a computer 
analysis for the prediction of the flow of hot gases, composed of a variety 
of exhaust products, through rocket and gas turbine components would 
support this program and also would be suitable for a wide variety of other 
applications. 

The Axisyrametric Diffuser Duct (ADD) code, developed by Anderson et 
al. (Ref, 1 and 2), has been shown to provide accurate and reliable viscous 
compressible flovpath calculations through small gas turbine components. 
It has been successfully applied to predict the flow in automotive gas 
turbine components (Ref. 2) such as that shown on Fig. 1, and swirling 
flow in a precombustion diffuser (Ref. 3). Vatsa et al«, (Ref. 4) has 
extended the analysis to the prediction of supersonic underexpanded 
cof loving jet flows and Anderson et al., (Ref. 3) has extended the analysis 
to treat supersonic internal flow with shock waves and expansion waves. 
Subsequent operational use of the ADD code at NASA Levis has identified a 
number of shortcomings in the code. The combined effects of low Reynolds 
number, large wall curvature, high swirl angles, and limitations in the 
turbulence models have taxed the capabilities of the code and produced 
stability and accuracy problems. The object of the present research 
program ie to: 1) refine both the algebraic and two equation turbulence 
models presently in the code to include tie effects of free stream 
turbulence and laminar/ turbulent transition, 2) provide accurate 
predictions of the thermodynamic and transport properties combustion 

product gases, and 3) improve the coordinate grid calculation for cases of 
long highly curved ducts. Following these code modifications, a matrix of 
test cases shall be calculated to: 1) identify numerical stability and 

computational accuracy problems occuring within the envelope of operation, 
2) correct problems where possible, and 3) establish grid selection 
criteria to reduce operational problems. 

DISCUSSION 

Both the algebraic and two equation turbulence models in tae ADD code, 
currently operational on the NASA Levis Cray computer, have been modified 
to include the effects of free stream turbulence and lamiaar/turbulent 
transition. The effects of free stream turbulence have been included in 
the algebraic turbulence model by modifying the outer layer Clauser 
constant using an empirical correlation developed by Yuhas (Ref. 6). 
Calibration tests of the accuracy of this model are shown on Figs. 2 and 3 
where the results are compared to the flat data of Blair (Ref. 7). Fig. 2 
shows a skin friction plot versus Reynolds number based on momentum 
thickness. As can be seen on this plot, the calculated skin friction for 
O.OZ turbulence follows very closely the skin friction law of Coles (Ref 
8} and is just slightly above the experimental data of Ref. 7. At the 3.0% 
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level the prediction is also slightly above the data> however the 
incremental change in skin friction has been predicted very accurately. 
Fig. 3 shows a comparison of the calculated velocity distributions with the 
experimental data of Ref. 7 for both the O.OZ and 5.0Z turbulence levels in 
universal coordinates. Also shovn on this plot is the law of the wall. 
Again» it can be seen that analysis shows excellent agreement with the 
experimental data in these calibration tests with the benchmark experiments 
of Blair (Ref. 7). 

The two equation turbulence model developed by Chien (Ref. 9), which 
is currently in the ADD code, is under evaluation. This model only 
requires a proper specification of the initial conditions for the turbulent 
kinetic energy and dissipation since the physics contained in the equations 
account for both free stream turbulence and natural transition. In the 
boundary layer, the initial conditions are specified in the same manner as 
described by Wang (Ref. 10). Outside the boundary layer, the turbulent 
kinetic energy (turbulent intensity) and dissipation (dissipation length 
scale) are specified. Although turbulent kinetic energy is frequently 
measured and therefore usually known, the dissipation is not. Therefore, a 
brief analysis was conducted to provide estimates of the dissipation length 
scales for different levels of turbulence. Two types of analyses were 
employed. The first method estimated the free stream dissipation from the 
rate of decay of turbulent kinetic energ> measured by Blair (Ref. 7). The 
second method estimated the dissipation from the turbulent kinetic energy 
spectrum, also measured by Blair, for different levels of free stream 
turbulence. A composite kinetic energy profile composed of both high 
frequency (dissipation length scale) and low frequency (integral length 
scale) correlations was fitted to the experimental data. Typical results 
of this latter procedure are sho%m on Fig. 4 which also shows the well 
known low frequency correlation of Von Karman. Based on these analyses, 
the dissipation length scale was found to be fairly censistant using both 
methods and was relatively constant for all levels of turbulence. 
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Figure 1 Computational Mesh for AGT101 Turbine 'niet Duct 
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FLOU-imUCfiD VIBRATION OF SSME MAIN INJECTOR LIQUID-OXYGEN POSTS 

S.S. Chen, J.A. Jendrzejczyk, ard M.W. Vambsganss 
Argomie National Laboratory 
Argonne, Illinois 60439 

The liquid-oxygen (LOX) posts are exposed to hot hydrogen flowing over the 
tubes on its way to the combustion chamber* Fatigue cracking of some LOX posts 
was observed after test firing of the SSMEs* A current design modification 
c^-^sists of attaching impingement shields to the LOX posts in the cuter row. 
a Ills modification improved the vibration/fatigue problem of the LOX posts, but 
resulted in an increased pressure drop that ultimately shortened the life 
expectancy of other components* 

A fundamental study of vioration of the LOX posts was initiated to 
understand the flow-induced vibration problem and to develop techniques to avoid 
detrimental vibrational effects with the overall objective of improving engine 
life. This paper summarizes this effort, including an assessment of the 
problem, scoping calculation and experiment, and a work plan for an integrated 
theoretical/experimental study of the problem* 

RESULTS OF WORK 

An assessment of the LOX post vibration problem has been performed to 
provide tae basis for an integrated theoretical/experimental study* Different 
excitation i!sechanisms were considered and the results are summarized as follows: 

1* Vortex Shedding: The synchronization of LOX post vib ation with vortex 
shedding is not possible in either drag or lift direction because of several 
reasons: relatively large mass-damping parameter, swirlers (spoilers) on the LOX 
posts, high Reynolds number, relatively small pitch-to-diameter ratio, and 
highly turbulent flow. 

2* Fluidelastic Instability: Based on a gap flow velocity oi about 

1670 ft/sec across the injector elements, the unshielded LOX posts are subjected 
to fluidelastic instability* Therefore, fluidelastic instability is identified 
as a potential mechanism that can contribute to LOX post damage in a short time* 

3* Inscability Due to Liquid Oxygen: The LOX posts will not become 

dynamically unstable due to the flow of liquid oxygen inside the LOX posts* 

4* Acoustoelastic Resonance: The lowest acoustic frequency is potentially 
close to the lowest natural frequency of LOX posts* However, the significance 
of acoustic resonance depends on acoustic excitation, which may be associated 
with vortex shedding or turbulence* Vortex shedding is unlikely to become an 
Important acoustic excitation* However, turbulence as well as other excitation 
sources have to be considered. 

5* Turbulent Buffeting: The detailed turbulence characteristics are not 
known at this time. Tui'‘alence burfetlng is a potential excitation mechanism. 

A scoping ( eperiment was designed and performed for a tube array simulating 
a typical grouping of LOX posts. The experiments were performed in a water 
loop* The tube arrangement is to simulate post 77 in row 13 and its surrounding 
tubes. Figure 1 shows typical tube accelerations as a function of flow 
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velocity* The movements of the tubes depend strongly on flow velocity. At low 
flow velocities, the motion is small. As the flow velocity increases, tube 
aceleration increases. When the flow velocity exceeds a certain value, tube 
response increases rapidly; this threshold flow velocity is called the critical 
flow velocity. The critical flow velocity obtained in the water loop was used 
to evaluate the critical flow velocity of the unshielded LOX posts. It is found 
that the gap flow velocity of 1670 ft/sec is in the unstable region# 

PERSPECTIVE/PLANS AND GOALS 

The dynamic response of the LOX post can be represented by the following 
equation 

[Mg + Mj]{Q} + [Cg + C^]{Q} + [Kg + Kj]{Q} = {G} , (1) 

where Q is a generalized displacement; and K. are the post mass, 

damping, and stiffness matrces, respectively; C^, and iC^ are the fluid added 
mass, fluid damping, and fluid stiffness matrices, respectively; G is the fluid 
excitation force. The key elements are C~, and G. 

At this time, there are no analytical or numerical methods available for 
calculating Cj and for a general array of cylinders in crossflow. An 
analytical solution in conjunction with experimental data is being explored to 
provide some additional insight into the characteristics of the fluid forces. 
Another method is the numerical technique based on the computational fluid 

dynamics (CFD) method. The advantages and disadvantages of different numerical 
methods are being reviewed. A computer code will be developed to calculate 
and Kj for the LOX post arrays* 

Fluid excitation forces will be measured in a test section for a tube array 
simulating the LOX posts. Fluctuating fluid force components are to be measured 
in the drag and^^Hft directions. Different grids will be used to simulate 
different upstream turbulence intensity. In addition, techniques will be 
developed to excite a tube in a tube array subjected to crossflow. At the same 
time the fluid force acting on the excited Lube as well as the surrounding tubes 

will be measured as a function of tube motion. Ftom these fluid forces, the 

required fluid force coefficients can be obtained. 

CONCLUDING REMARKS 

The LOX post array subjected to hot gas flow is a very complicated 
system. The detailed flow field is not known at this time, so a detailed 
evaluation of the LOX post vibration potential is not possible. The cential 
issue in predicting the LOX post response is the fluid excitation force and 
motion^ependent fluid force. This remains a difficult issue. The work 
initiated by Lewis to develop techniques and to conduct tests to obtain the 
necessary fluid force is essential in the understanding of the vibration of the 
SSME LOX posts and in the development of techniaues to avoid detrimental flow- 
induced vibration. 
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OVSRVIEU OF STRUCTURAL RESPONSE: PROBABILISTIC STRUCTITRAL ANALYSIS 


Christos C. Chamic 

National Aeronau'^ics and Space Administration 
Lewis Research Center 
Cleveland » Ohio 44135 


Advanced analysis methods are required to predic!: accurately the struc- 
tural response (static, transient, cyclic, etc.) and acconq>anying local 
stresses in space propulsion system components operating in a fatigue environ- 
ment consisting of complex thermal and mechanical load spectra. Advanced 
three-dimensional analysis methods must be formulated to accommodate the 
limited deterministic and probabilistic complex loading spectra as well as th^ 
anisotropic, thermoviscoplastic behavior of the material. Presently three- 
dimensional, deterministic, anisotropic, thermoviscoplastic analysis methods 
are not available and recourse is made to incremental, iterative analysis 
capabilities available in geaeral-purpose structural analysis programs. Rec- 
ognizing the inq)ortance and necessity of Ceveloping solutions for con^lex 
three-dimensional inelastic problems associated with aeronautical propulsion 
system components, Lewis is conducting intensive research under the HOST (Hot 
Section Technology) and Research and Technology (R&T) Base progiams to 
develop the requisite deterministic analysis methodology for the solution of 
these problems. The deterministic, three-dimensional, inelastic analysis 
methodology developed under these HOST and RAT Base programs will be aug- 
mented to accommodate the complex probabilistic loading spectra, the thermo- 
viscoplastic material behavior, and the material degradation associated with 
the corrosive and oxidative envi^^onments of space propulsion system structural 
components representative of the space shuttle main engine (.^SME) such as 
turbine blades, transfer ducts, and liquid-oxygen posts. 

The probabilistic approach to structural response consists of the follow- 
ing program elements: (1) composite load spectra, (2) probabilistic structural 

analysis methods development, (3) probabilistic finite element theory, 

(4) probabilistic structural analysis application, (5) stnictural tailoring of 
turbopump blades, (6) unified theory of dynamic creep buckling/ ratcheting, (7) 
creep buckling/ratchecing analyzer, and (8) nonlinear COBSTRAN development. 
Research activities on all of these program elements (except the creep 
buckling/ratcheting analyzer) are under way. Program participants, program 
element description, and progress are the focus cf this session. 
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SSME STRUCTURAL DURABILITY 


STRUCTURAL RESPONSE 


BACKGROUND: 

0 LOADING CONDITIONS ARE NOT DETERMINISTIC IN GENERAL 

0 MATERIAL THERMOMECHANICAL BEHAVIOR DIFFICJVT TO CHARACTERIZE IN SSME SERVICE 
ENVIRONMENTS 

0 CONTINUING EVOLU'ON OF NEW MATERIAL^ FOR USE IN SSME 

0 CONTINUING INCREASES IN PERFORMANx^t REQUIREMENTS OF SSME 


SSME STRUCTURAL DURABILITY 

STRUCTURAL RESPONSE 


RESULT: 

0 DETERMINISTIC ftTHODS ARE INADEQUATE TO DESCRIBE THE STRUCTURAL RESPONSE OF SSME 
CRITICAL COMPONENTS WHERE; 

0 THE LOADING CONDITIONS ARE M DETERMINISTIC 
0 THE MATERIAL THERMOMECHANICAL BEHAVIOR IS NOT FULLY CHARACTERIZED 
0 THE DEMANDS FOR NEW MATERIAL AND INCREASED PERFORMANCE IS EVaUTIONARY 
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SSME STRUCTURAL DURABILITY 


structural response 


NEED/OBjECTiVEt 


0 DEVELOPMENT OF PROBABILISTIC METHODS 


0 LOADING CONDITIONS DEFINIflQN 


Q material behavior DFSCRIPTION 


0 STRUCTURAL ANALVSib PREDICTION 


FOR CRITICAL SSME STRUCTURAL COMPONENTS IN ATTENDANT SERVICE ENVIRONMENTS 


SbME SIRUOTURAL UUHAUlLllY SIRUCIUNES RESEARCH AtllVlTIES 


PROGRAM ELEMENT /ACTIVITY 



GOMPUSllE LOAD SPECIRA 


PROOHA/i SCHEDULE 

‘ f — — ^ 





PROOAUILISIIC SIRUCIURAL 
ANALYSIS ME I MUDS 


V.AHIAdONAL AITROACH lU GRANT 

P^iOnAblLISTlC FINIIE ELEMENIS 


PROBABILISTIC SIRUCIURAL 
A/ML YS IS 


IN-liOUSE 


S TRUC j URAL ANAL YS IS / 1 AILOKING IN-HOUSE 
OF iURBOPUHP BLADES 


UNIFJEU THEORY OF DYNAMIC 
UiEEP OUCKLING/RAICHEUNG 


CREEP BUCXLIN&/RA1CHEUNG 
ANALYZER 


NONLINEAR COBSTRAN 


IN-HOUSE 
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S)RUC1UR£S RESEARCH ACimUES ON SSHE S1RUCTURAL UURABILIIY 


RtSCARCH ACTIVITY TYPE 

COMPOSITE LOAD SPECTRA CONTRACT 

PKOBABUISTIC STRUCTURAL CONTRACT 

ANALYSIS NETHODS 

VARIATIONAL APPROACi! TO 6RANT 

PROBABILISTIC flNITC ELEMENTS 

PROBABILISTIC STRUCTURAL IN-HOUSE 

analysis 

STRUCTURAL ANAL YSIS/TAIL0RIN6 IN-HOUSE 
Of TURBOPUHP BLADES 


fl&iEf MSCRIPTIQH 


PKOBABLISTIC METHODS TO DEFINE THE INDIVIDUAL 
AND COMBINED LOADS ON SSME COMPONENTS 


PROBABILISTIC DETERMINATION OF GLOBAL/LOCAL 
NONLINEAR DISPLACEMENTS/STRESSES/ STRAINS 


APPLICATION OF EXISTING PROBABILISTIC METHODS 
TO STRUCTURAL ANALYSIS OF SSME TYPE COMPONENTS 


ADArlATlON Of STAEBL TO SSME TURBOPUMP BLADES 


DEVELOPMENT OF ADVANCED PROBABILISTIC 
STRUCTURAL MECHANICS METHODS 


UNIFIED THEORY OF DYNAMIC GRANT 

CREEP BUCKLING/RATCHETING 


DERIVATION AND INTEGRATION OF MECHANISTIC 
MODELS OF VARIABLES INFLUENCING CREEP 
BUCKL1N6/RATCHETING 


CREEP BUCKLING/RATCHETING CONTRACT 

ANALYZER 


DEVELOPMENT OF A COMPUTER CODE TAILORED FOR 
CREEP BUCKLING/NATCHETING 


NONLINEAR C06STRAN 


IN-HOUSE ADAPTATION OF COBSTRAN TO ACCOMMODATE THE THERMO' 
VISCOPLASTIC BEHAVIOR OF SSME TURBOPUMP BLADES 


STRUtlURAL REbPONSE OENEFilS 


0 SIKjRl TERM: 

0 HORE Ati^iRATE TIME-PHASING OF LOADS 
0 IMPROVED ACCURACY OF PHEDICUD SIRESS/SIRAIN 

0 WllH ASSESSED PR00A81L11Y AND CONFIDENCE LEVEL 
0 SPLtlFlC GUIDELINES FOR DESIGN MOOIFICAIIONS 10 INCREASE SERVICE LIFE 

0 LONG TERM: 

0 PREDICIIVE METHODS 10 DESIGN SSME CRIIICAL COMPONENTS WITH SPECIFIED 
0 SIRUCIURAL DURABILITY 

0 SERVICE LIFE 

0 PROBABILITY 

0 CONFIDENCE LEVEL 
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COMPOSITE LOADS SPECTRA FOR SELECT SPACE PROPULSION STRUCTURAL COMPONENTS* 

J.F. Newell 

Rockwell International 
Canoga Park, California 91304 

Rocket engine technology continues to demand higher performance with 
lighter weight components that have man-rated reliability requirements. These 
requirements have yielded higher operating pressures, temperatures, and tran- 
sient effects as well as markedly increased mechanical vibration and flow- 
related loads. The difficulty in installation, cost, and potent* 1 for 
destroying an engine severely limited the required instrumentation and measure- 
ments to adequately define loads of key components such as turbine blades. 

Also, accurate analytical methodologies for defining internal flow-related 
loads are just emerging for problems typically found in rocket engines. The 
difficulty of obtaining measured data and verified analysis methodologies has 
led to the probabilistic load definition approach of this contract. Rocketdyne 
has the state-of-the-art methods for defining loads and analysing structural 
problems and has teamed with Battelle to utilize their expertise in the devel- 
opment of probabilistic models. 

The objective of this program is to develop generic load models, with 
multiple levels of progressive sophistication to simulate the coioposite (com- 
bined) load spectra that are induced in space propulsion system components, 
representative of space shuttle main engines (SSME), such as transfer ducts, 
turbine blades, and liquid oxygen poses. These models will be developed 
through describing individual loading conditions using state-of-the-art prob- 
abilistic methods combined with deterministic coefficients related to key 
engine variables. The combinations of these loading conditions will use a 
similar approach to syntnesize the composite load spectra. The use of key 
variables, such as power level, flow, pump speed, or weight, allows for generic 
load definition for scaling and prediction of component loads. 

The second approach in the option portion of the contract will consist of 
developing coupled models for composite load spectra simulation which combine 
the (mechanistic) models for thermal, dynamic, acoustic, high-pressure, and 
high rotational speed, etc. Probabilistic synthesis using statistically vary- 
ing coefficients will be used to obtain the composite load spectra. 

The unified theory required to combine the various individual load simu- 
lation models (hot-gas, dynamic, vibrations, instantaneous-position, centri- 
fugal, etc.) into composite load spectra simulation models will be developed 
under this program. Results obtained from these models will be compared with 
available numerical results, with the loads induced by the individual load 
simulation models, and with available structural analysis results from inde- 
pendent analyses and tests. These theories developed under both approaches 
will be further validated with respect to level of sophistication and relative 
to predictive reliability and attendant level of confidence. 


*Work performed under NASA contract NAS3-24382. 
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A computer code configured as an expert system will be developed that 
incorporates the various individual and composite load spectra models to con- 
struct the specific load model desired. The approach is to develop and deliver 
the computer code at four intervals in the contract. The first delivery will 
be an initial code for turbine blade loading. Each subsequent delivery will 
add sophistication to the component probabilistic load definition and the 
decision-making processes as well as installing a new set of loads for an 
additional component (e.g., liquid oxygen post). This allows for on-going 
evaluation and usage of the system by both Rocketdyne and NASA starting with 
the first code delivery. 

The initial efforts for this program are proceeding with the probabilistic 
model theory. The surveys of individual loadings, combined loadings, proba- 
bilistic methods, and mission history profiles are in progress. The liquid 
oxygen post and turbine blade survey efforts have been completed. The survey 
of appropriate probabilistic methodologies applicable to the rocket engine has 
proceeded to testing and evaluating methodologies with typical SSME test data. 
This effort is summarized in a separate presentation. 

The simulation of individual load models for turbine blades has been 
started so that the initial turbine blade code can be delivered by October. 

This includes both relating measured test data to the individual turbine blade 
load coinponent-like pressure, temperature, or centrifugal loads and developing 
probabilistic coefficients from engine data. 

The code development effort has just been starced in parallel with the 
other tasks. Defining the overall approach and knowledge base requirements is 
the first task. 
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Rockwell 

International 

R9cl(»ldyn« Dhrfskm 




700 603 


GENERIC LOADS 


• COf.SIDER ALL POTENTIALl.V SIGNIFICANT LOADS 

• SSME SUBSET PRIMARILY USED 

• LOADS APPLICABLE TO LO2/LH2 ENGINES 

• LOADS RELATED TO COEFFICIENTS AND KEY PARAMETERS 

• COMBINATION OF DETERMINISTIC AND PROBABALISTIC VARIACLES 

• KEY >•' ^AMETERS AFFECTING LOADS INCLUDE 

• ENGINE CONFIGURATION • FLOW INTERRUPTIONS 

• THRUST • PUMP TORQUE 

• i'.cIGHT • »*UMP SPEEDS 

• INLET CONDITIONS • ETC 

• COMPONENTS 

• LOX POST • TRANSFER DUCT 

• TURBINE BLADE • OTHER - MCC LINER, HPOTP DUCT, NOZZLE FEEDLINES 
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DATA ANALYSIS 


8ASIC EQUATIONS 
AND CCEFICIENTS 
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-• 

TRANSIENT 

ANALYSIS 
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• DETERMINISTIC 
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CORRE'^TION 
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COMPONENT 
RELATIONSHIPS 
(TURBINE BLADE) 


HPFTP TURBINE OPERATING CONDITIONS 
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FLOW INTERRUPTIONS 
A 13 struts 
9 41 NOZZLES 
C 63 BLADES 
O 39 VANES 
E 59 BLADES 
F 12 VANES 






3885 PSIA 
10OO®R 



LINER 

4943 PSIA 
4950 PSIA 




Hocliwall 
International 
RocHtMiFt'* DK1*I0( 


ORiG«flAL Pf-O 
OF POOR C J.. 









«*S,PS»A,OEQ R,RPM. POUNDS 


HIGH PRESSURE FUEL TURBINE qf pc o T ri 
TRANSIENT CHARACTERISTICS 



TEMPERATURE OVNAMIC LOAD I TEMPlIRATURF 

CRITICAL CRITICAL CRITICAL 



TURBINE BLADE LOADS 



Rock^l 

In.efnetHjoal 

AecfcstdpM DlrMdn 


• AIRFOIL 

• centrifugal 

• PRESSURE 

• BLADE TEMPERATURE 
OVNA .-’C 4P 

• IMPiNGrMENT LOADS 

• OAMPER/FRICTION 

• TIP RUBBING 

• SHANK 

• PRESSURE 

• METAL TEMPERATURES 

• EXPERT OPINION INPUT REQUIRED 

• LACK OF MEASURED DATA 


KEY VARIABLES 
S"-ED 

TURBINE ANALYSIS 
TURBINE/THERMAL ANALYSIS 
SPEED. INTERRUPTIONS 
UPSTREAM DEBRIS POTENTIAL 
GEOMETRY, WHIRLIGIG TEST 
DESIGN 


COOLANT FLOW AND THERMAL 
ANALYSIS 
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489^773 


LOAD SPECTRA CODE 


• iNlTIAL PLANNING STAGE 

• OERNING KNOWLEDGE BASE REQUIREMENTS 

• INDIVIDUAL LOAD PROBABILISTIC FUNCTIONS 

• COMPOSITE LOAD SIMULATION TECHNIQUES 

• PROBABIUSTfC METHODS AND SOLUTION ALGORITHMS 

• SSME BASELINE INFORMATION 

• KEY VARIABLES MISSION HISTORY PROFILES - GROUND 
TEST AND FLIGHT 

• ENGINE SYSTEM PACKAGING AND SIZING 

• ENGINE CYCLE AND COOLANT TECHNIQUES 

• COMPONENT CHARACTERISTICS 

• MEBl OVERALL GOALS OF CONTRACT 

• CODE SOPHISTICATION ENHANCED IN EACH OEUVERED 
MODULE 

0 

1 


INTERRELATIONSHtPS OF CONTRACT EFFORTS 
SSME AND PROBABILISTIC UTILIZATION 


• MATURE PROBABILISTIC METHODOLOGY INTEGRATED 
AND USED aN SSME AND FUTURE ENGiNES 


• PROBABILISTIC STRUCTURAL 
ANALYSIS CONTRACT 


• LOADS SPECTRA DEFINITION AND MODEL 
FOR PREDICTION AND ASSESSMENT OF 
STRUCTURAL DURABILITY OF SPACE 
PROPULSIOF’ COMPONENTS 


• CURRENT SSME DEVELOPMENT WORK 
% dynamic DATABASE 
% DECISION ANALYSIS OF TEST DATA 
« BASED ON MEASUREMENTS AND 
HARDWARE CONDITIONS 
% BOTh ROCKETOYNE AND MSFC EFFORT 


• CURRENT OR NEW 

NON PROBABILISTIC ANALYSIS 
METHODOLOGIES 


• SSME OPERATION 

• FABRICATION 

• TESTING 

• OPERATIONAL FLIGHTS 

• PRODUCT IMPROVEMENTS 

• PHASE II 

• MSFC TEST BED 

• ETC 


Rockwell 

Intematfof»al 

Rochtldyn* Division 


75 







N85 -27954 

COMPOSITE LOADS SPECTRA FOR SELECT SPACE PRCPULSIOR STRUCTURAL COMPONENTS - 

PROBABILISTIC LOAD MODEL DEVELOPMENT* 

Robert Kurth 

Battelle Columbus Laboratory 
Columbus , Ohio 43201 

This effort is in support of the development of an expert system of com- 
puter codes to predict the loads on select structural components of a space 
propulsion engine. The development will be based primarily on the space 
shuttle main engine (SSME) test data base. Because of random variations of the 
many different sources of the loadings on the selected structural components 
and transients, a probabilistic approach to the problems has been adopted. The 
goal of this task is to characterize all of the individual sources of loads at 
critical structural locations, such as the turbine blades, the transfer ducts, 
and liquid oxygen posts, using state-of-the-art probabilistic methods with 
varying levels of sophistication. 

The second phase of this work is the development of a composite load model 
based on a probabilistic synthesis of the individual load model previously 
developed. This model will be based on the stochastic combination of the load 
variables and not on the physical process for the combination of the individual 
loads to the composite load seen by the selected structural components. 

ne probabilistic load model effort has been started with an initial sur- 
vey c f available methods and an assessment of these methods using typical SSME 
data. The results of the first year of study will produce a model applicable 
to turbine blades in an engine represented by the SSME. The survey of appro- 
priate statistical methods has been completed for this phase of the work and 
the assessment of these methods is in progress. Th^ engine measurements 
related to i-urbine blade loads are under study primarily to support the first 
year code development effort. 




page blank not Fn.MEO 


^Work performed under NASA contract NAS3-24382. 
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SURVEY OF AVAILABLE METHODS 


• LllERATURt SURVEY CONDULUD 

• MANY INVESTIGATIONS HADE BUT NUHBER OF DUFEREN" METHODS LIMITED 

• MOST IMPORTANT MEIHODS INCLUDE 

•ANALYTIC 

•FIRST ORDER SECOND MOMENT 
•BARRIER CROSSING (FIRST PASSAGE! 

•DISCRETE SIMULATION 






OBaneiie 


ANALYTIC METHODS 


• MATHEMATICALLY COMPLEX 

• ONLf WORK FOR IIMHID SET OF DISTRIBUTIONS 

• ' IHITED APPLICATION TO SSHE 




Ror^wvll 


OBaifeiie 
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FIRST ORDER SECOND MOMENT METHODS 


• IN USE PRIMARILY IN SlRUCTURAL ANALYSIS 

•RACKWITZ-FIESSLER 

•CHEN-LIND 

•WIRSCHING 

• PROVIDE THE PROBABILITY OF A LIMIT STATE BEING EXCEEDED 
•REQUIRES MANY APPLICATIONS 10 DEFINE PROBABILISTIC NATURE OF LOADS 

• coordinate 1RAN$^0RMAT10NS AND APPROXIMATIONS NEEDED TO HANDLE 
NON-NORMALLY DISTRIBUTED VARIABLES 

• LIMITED USE FOR INDIVIDUAL LOADS BUT HAY BE U;tFUL FOR COMPOSITE LOAD MODEL 

OBaiteiie 


DISCRETE SIMULATION METHODS 


• WILL BE NEEDED fOR THE EXPERT CODE DEVELOPMENT AT LEAST TO VERIFY THE 
CALCULATIONS 

• MOST IMPORTANT METHODS ARE 

•MONTE CARLO 
•MARKOV CHAIN 

• DISCRETE PROBABllITY DISTRIBUT'ON 

• rascal 




noct-«w^ 


OBatfeiie 
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OWQWAL BARRIER CROSSING METHODS 

OF POOR Q^’ALITY 


• r'ROVlDf.S hN FypFCILD NUMBCR. OR RATt OF. CROSSINGS OF A GIVEN LEVEL, I E . 
THE BARRIER 


• CAN BF U‘^'-D TO DFTF.RHlNf THE PR03ABII «TY OF ANY GIVEN IFVFl BEING EXCEEDED 


• RELATIVELY STRAIGHIFORWARD Ij APPLY 




RocH*v<( 


OBanei'e 


CHARACTERIZATION OF CURRENT LOADS 


• OESlRf P PFSULTS ARF 

• ''SUPATF OF THE MEAN 
•ESTIMATE OF CREDIBLE VARIATIONS 

• MOST SUITABLE PDF EORM 

• stationarhy test 

•CORRELATION EFFECTS. IF ANY 

• FREQUENCY OF ( OADS 

•DISTRIBUTIONAL TAMILY IS CHOSEN BY DISTRIBUTION SHAPE. I.E . 
SKfWN. SS AND KURTOSIS 

• OISTRfBf/TION PARAMEIIRS OBTAINET B' MOMENT MAT"* ’SS 
Chi SQUARE TESl US^O TO CMFCK 




nocV«««N 

IntumaHoml 


OBdiieiie 
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PROB/>SIUSTIC FUNCTION DEFINITION 


• GOAl IS TO DEMNf IMf I OAD LEVEL, Lj. AND ITS PROBABILITY OE OCCURRENCE. 


• THE LOAD LEVII. Lj IS DEMNED BY A FUNCTION WHICH IS A COMBINATION OF 
DETERMINISTIC AND PROBABILISTIC METHODS 

Li(T) Fo (Xi. X?. .X(,J 0Fp(Yi. Y'^ Yh) 

Xi’ DtTFRMINISTIC INPUTS 
Yl PROBABILISTIC INPUTS 
Fd DETEfIMINIsnC FUNCTION 
PROBABILISTIC FUNCTION 


• SOME. OR AIL, OF IHF Xj'S AND Y^'S MAY BE THI SAME 


*» 


Rocf«M« 

Intentional 

o m n n 
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PR03A3ILISTIC FUNCTION DEFINITION 


rt THE FUNCTIONS AND fp ARE NOT NECESSARILY ANALYTIC, E.G. THEY HAY 
REPRESENT COMPUTER CODES 


• THE COEFFICIENTS AND lOMPILXITy OE THE FUNCTION I . (T) WILL BE 
ETERMINED BY THE ACCURACY REQUIRED BY THE ANALYST 


• TXPERT OPl'^ION 15 CONTAINED IN THESL ^'"^OLATIONS 


411^ Mocl-***! 

Wen«llon*l 


OBaneiie 
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PROBABILISTIC FUNCTION DEFINITION 


• THE PROBABILIIY OF iHF LOAD L^U) OCCURRING IS GIVEN BY 

Paj(T)} = PCY^lPlYj) ..P(Yj,) 

• IVIS PROBABILITY CAN BE ONE 

• IF Y , 1$ DtPLHDtNT ON THEN U'NDUIONAL PROBAB\l UY IS OSLO 

J “ 

• TO ADDRESS THE VARIATTO: ABOUT THE LOAD Uvfl I^IT) BARRIER CROSSING 
METHODS ART USFD 



OBa!felie 


PROBABILISTIC FUNCTION OUTPUT 


• the PROBABILISTIC FUNCTION OUTPUT IS AN ARRAY OF DISCRETE LOAD LEVELS WITH T.TEIR 
ASSOCIATED PROBABILITY OF OCCu^StNCE 

R = t(L^. P^), (L^r P 2 ). .. (1^. P,|)l 


• THE reliability (ACCURACY) OF THE OUTPUT CAN BE CONTROLLED BY INCREASING OR 
DECREASING N 


• THE VECTOR or ORDERED PAIRS R IS A DISCRETE RTPREsENTA '-.N Of THE DENSITY 
FUNCTION FOR THE I OAD AND CAN BE USED TO CAI CULATI MEANS, VARIATIONS, AND SO FORTH 



OBditelie 
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Or \ 


STATUS OF PROBABILISTIC LOAD MODEL 
DEVELOPMENT 


• sumt COMPLUID 

• STAIISTICAL ANALYSIS PACKAGE. COHPLETED 

• BARRILR CROSSING ANAIYSIS MEiHOD COMPUTED 

• TURBINE BLADE LOAD ANALYSIS STARTED 
•JOINT DENSITY FUNCTIONS CALCULATED 

• PROP.AH’i ISMC fUNCTlON DEVELOPMENT 3EGUN A«t TESTED FOR 65t AND 
IC'U R ' IVU S 




OBaneiie 
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N85-27955 

PROBABILISTIC STRUCTURAL ANALYSIS THBORY DBVBLOPMBNT* 

O.H. Burnsic'e 

Souttmest Research Institute 
San Antonio, Texas 7828A 

The objective of the PSAM project is to develop analysis techniques and 
computer programs for predicting the probabilistic response of critical struc- 
tural components for current and future space propulsion systems- This tech- 
nology will play a central role in establishing increased system performance 
and durability. 

Southwest Research Institute has assembled a project team from industry 
and the university community whose technical expertise and experience make 
them uniquely qualified for the PSAM effort. The participants on the team 
are: Sou^:hwest Research Institute - project management and solution strategies 

develC}' V MARC Analysis Research Corporation - code development for probabi- 
listic unite element methods; Rocketdyne Division, Rockwell International 
worp. - Space Shuttle Main Engine (SSME) design and hardware experience; Prof. 

Paul Wirsching (University of Arizona) - probabilistic and reliability methods; 

Prof. Gautam Dasgupta (Columbia University) - stochastic finite elements; 
and Prof. Satya Atluri (Georgia Tech) - advanced material constitutive models. 

The PSAM program was initiated in October 1984 and is a four-year effort. 

The first year's technical activity is concentrating on probabilistic finite 
element formulation strategy and code development. Work is also in progress 
at Rocketdyne to survey critical materials and SSME components such as turbine 
blades, liquid oxygen posts, and transfer ducts. This survey will provide 
information on the operating environment and failure modes encountered by 
space propulsion system components, as well as verification cases for the 
computer codes developed in PSAM. 

The probabilistic finite element computer program NESSUS (Numerical Evalua- 
tion of Stochastic Structures Under Stress) is being developed by MARC Analysis 
Research. MARC's code development activities for NESSUS are presented elsewhere 
in these proceedings. The probabilistic FEM code development is scheduled 
for three years. The final probabilistic code will have, in the general case, 
the capability of performing nonlinear dynamic analysis of stochastic struc- 
tures. Stochasticity can be found in the loads and in the structure due to 
material properties, geometry, and boundary conditions. The goal of the first- 
year effort is to develop a probabilistic FEM code for NASA-LeRC which has 
the capability to analyze linear elastic structures under static and dynamic 
loading. With this code, the dynamic response of components such as turbine 
blades can be computed. Nonlinear enhancements will be made in the final 
two years of the finite element task. 

The remaining two technical tasks. Approximate Methods and Advanced Methods, 
will begin in the second year of the PSAM program. It is the goal of the 
Approximate Methods effort to increase problem solving efficiency relative 
to finite element methods by using energy methods to generate trial solutions 
which satisfy the structural boundary conditions. These approximate methods 
will be less computer intensive relative to the finite element approach because 


*Work performed under NASA contract NAS3^24389. 

85 

Ai VjLXSii NOT ITLMFD 





of the reduced number of degrees of freedom and also the reduced number of 
random variables. The object of the Advanced Methods task is to stek new 
algorithms which can be applied to probabilistic structural analysis, rather 
than to refine existing ones. While selection of such methods has not been 
made, several candidates are under consideration. These Include using weighted 
residual techniques for deriving the structural equations and the boundary 
element method applied in a probabilistic format, together with artificial 
intelligence/expert systems to aid the engineer in performing probabilistic 
analyses. The Advanced Methods task will be conducted over the final three 
years of the PSAM program. 

Three approaches have been identified for computing the probabilistic 
structural response. These are fast probability integration (FPI) methods 
used in reliability analysis, methods of moments, and Monte Carlo simulation. 
Both FPI and moment methods are computationally efficient, but do involve 
certain approximations. Monte Carlo gives "exact** results, but requires a 
large number of realizations to obtain high accuracy. For large systems with 
many random variables, the Monte Carlo method does not appear practical because 
it is so computer intensive, especially if information about the tails of 
the distributions is desired. However, it is a valuable tool when applied 
to smaller systems for judging the accuracy of other methods. 

In summary, the analysis methods and computer codes developed under PSAM 
are expected to provide the tools for both enhancing the reliability of space 
propulsion system designs and advancing the state of probabilistic structural 
mechanics. 
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ORIGINAL PAG2 IS 
OF POOR QUALITY 


PROfiABIclSTIC STRUCTURAL ANALYSIS NETHOOS (PSAH) 
FOR SELECT SPACE PROPULSION SYSTEW COMPONENTS 


NASA-LeRC Controct Ik.. NAS5-2<i389 
Si«l Project No. 06-8T39 


- Progr'* Objective - 


Develop Analysis Tectmlques and Conputer 
ProgroRS for Predicting the Probabilistic 
Response of Crltlcol Structural Conponents 
for Current and Future Space Propulsion 
SysteoB 


PSAH PROJECT TEAM 


Priie CoaLfOCLor 

0 Southwest Research Institute 
Son Antonio, Texos 


Subcoatroctors 

0 NARC Ano lysis Research Corporation 
Polo Alto, Colifornlo 

0 Rocketdyne Division, Rockwell Intemotlonol Corp. 
Cnnogo Pork, Colifornlo 

0 University of Arlzono <Prof. Paul Nlrsching) 

Tucson, Arizona 

0 Georqlo Institute of Technology (Prof. Sotya Ailurl) 
Atlontn. Georgio 

0 Professor Goutcw Dosgiioto 
ColcJRblo University 
New York, New York 


87 



PROBABILISTIC STRUCTURAL ANALYSIS METHODS 
PROGRAM SCHEDULE 


TASK DESCRIPTION 


1 FINITE ELEMENT METHODS 


2 APPROXIMATE methods 


3 ADVANCED METHODS 


4 REPORTING 


i=Y 86 EV 87 FY 88 


5, MANAGEMENT 



Or POOH O'-"'-*' 


STOCHASTIC 

MAftERlAL 

PROPERTIES 




k STRESS 
\ CONCENTRATiaN 


STOCHASTIC 

THERMOMECHANICAL 

LOADS 


UNCERTAIN 

BOUNDARY 

CONDITIONS 



DYNAMIC 

RESPONSE 


uncertain 

STRUCTURAL 

SHAPE 


THE PROBABILISTIC 
STRUCTURAL ANALYSIS PROBLEM 
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PROBABILISTIC LOADING CONDITIONS FOR PSAH 


r. Pressure 
0 TDennol 
0 Centrlfugol 
0 Periodic 
o Acoustic 
0 Impulsive 


STRUCTURAL NONL ! NEAR 1 TIES IN PSAN 


0 Large Strolns 

0 Large Displacements and Rotations 

0 Natertal Behavior 

(Plosticity and Thermovlscoplostlclty) 
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CONCEPT OF PROBABILISTIC 
STRUCTURAL ANALYSIS 



INPUT PDF 

(Load, Material. Geometry. Boundary Conditions] 




STRUCTURAL MODEL 


o 



OUTPUT PDF 

(Displacement, Strain. Stress] 


osioiWAi. 

OF POOR QU 
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.•t. 


c,. 


1 • STRUCTURAL PARAMETERS ARE RANDOM 
BUT SPATIALLY HOMOGENEOUS 



THREE LEVELS OF SOPHISTICATION 
FOR PROBABILISTIC FINITE ELEMENT ANALYSIS 


STRATEGIES FOR COIVUTING THE 
PROBABILISTIC STRUCTURAL RESPONSE 


0 Post Probability Integrotlon Technlcjes 
0 Wethod of Homenis 
0 Nonte Corio Simulation 
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PSAM INTEGRATES EXISTING TECHNOLOGIES 
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N85-27956 

PROBABILISTIC FINITE ELEMENT DEVELOPMENT* 


Jop Nagtegaal 

MARC Analysis Research Corporation 
Palo Alto, California 94306 


The probabilistic finite element computer program NESSUS (Numerical 
Evaluation of Stochastic Structures Under Stress) is being developed by the 
MARC Analysis Research Corp. for the analysis of critical structural components 
for reusable space propulsion systems. This program has evolved from the HOST 
code, developed by MARC for Pratt & Whitney Aircraft Co. under NASA contract 
NAS3-23697, and which is now in use for the analysis of hot-section turbine 
components at Pratt & Whitney and the NASA Lewis Research Center. 

The development of the NESSUS code is scheduled to take 3 years. First 
year efforts involve the formulation of the probabilistic analysis strategy and 
the development of a probabilistic linear analysis code. The ultimate goal of 
the 3-year program is the development of a finite element code capable of per- 
forming nonlinear dynamic analysis of structures having stochastic material 
properties, geometry, and boundary conditions and subjected to random loading. 

Three levels of sophistication are envisioned for the stochastic descrip- 
tion of the structural problem, namely: 

Level 1: Homogeneous random variable for stiffness, mass, damping, and 

external loading 

Level 2: Stochastic characterization of variables at the element level, 

with specified interelement correlations 

Level 3: Stochastic interpolation of variables within a finite element 

Two alternative probabilistic analysis methods will be developed, allowing 
for all three levels of modeling sophistication: 

(1) A probability integration method, providing a direct estimate of the 
reliability of the structural configuration under study 

(2) A simulation method, providing the means to verify the results 
obtained with method 1 

Part of the formulation developed by Southwest Research Institute for the 
level 1 linear problem is currently available in a form that can be readily 
implemented in the NESSUS code. SwRI is continuing work on the formulation for 
the level 1 dynamic analysis. Which should be available in the near future. 
Progress is being made in the formulation for the level 2 static prol lem. The 
major issue at this level remains the choice of an efficient algorithm for the 
evaluation of the derivatives of the stiffness, mass, and damping matrices with 
respect to each random variable present in the analysis. An approach to bypass 
this problem by efficiently extracting the static response of a perturbed 


*Work performed under NASA contract NAS3-24389. 
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system without the need for explicit differentiation has been proposed by MARC. 
This method is based on a modified Newton-Raphson iteration scheme and could 
make use of some of the special numerical algorithms developed for the HOST 
program. MARC is presently testing this approach on a realistic problem while 
efforts are continued in the development of a similar formulation for the per- 
turbed eigenvalue problem. The extension of this method into the nonlinear 
regime appears promising, especially in systems offering parallel computation 
capabilities. 


NESSUS 

liUMBRlCAL EVALUATION OF 
STOCHASTIC STRUCTURES UNDER STRESS 


0 A COMPUTER CODE FOR THE PROBABILISTIC FINITE ELEMENT ANALYSIS 
OF CRITICAL COMPONENTS FOR SPACE PROPULSION SYSTEMS 


0 EVOLVED FROM THE QQSI CODE DEVELOPED UNDER NASA CONTRACT 
NAS3-23697 FOR THE ANALYSIS OF HOT-SECTION TURBINE PARTS 


0 THREE YEAR DEVaOPMENT PROGRAM. RESULTING IN A CODE CAPABLE 
OF PERFORMING NONLINEAR STATIC AND DYNAMIC ANALYSIS OF A 
STRUCTURE WITH STCiCHASTIC MATERIALS. GEOMETRY, BOUNDARY 
CONDITI(»«iS AND LOADS. 
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NESSUS 


FINITE ELEMENT LIBRARY 




3D BRICK 



3D SHELL 



LINEAR BEAM 


0 ALL ELEMENTS INTEGRATED NUMERICALLY 
0 LINEAR LAGRANCIAN INTERPOLATING FUNCTIONS 
0 SELECTIVELY REDUCED INTEGRATION 


NESSUS 

NONLINEAR ANALYSIS FEATURES 


0 NONLINEAR PROBLEM MAY BE SOLVED WITH 

(a) PURELY ITERATIVE METHODS 

(b) INCREMENTAL-ITERATIVE METHODS 

o THREE TYPES OF INaASTIC CONSTITUTIVE MODELS 

(1 ) SECANT aASTICITY SimpUTied PUaUdty 

(u ) VON MISES PLASTICITY Radial Raturo. AasociiUve Flow 

(U) CRUP-PLASTICin Walker s Model 


0 USER SUBROUTINES FOR TEMPERATURE-DEPENDENT ’ i aTERI AL 
PROPERTIES. ANISOTROPY. CREEP. W0RIC-HAR1»F':NG. ETC 
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NESSUS 


DYNAMIC ANALYSIS PEATUR£S 
for damped or uodamped systems 


1. EIGENVALUE EXTRACTION 

0 SUBSPACB ITERATION 
0 POWER SHIFTS 

2. LINEAR DYNAMICS 

0 MODE SUPERPOSITION 

(i ) FREE VIBRATION 

(U ) TRANSIENT forcing FUNCTION • 

(iii) HARM(Y. NODAL LOADING * 

(iv) HARMONic BASE EXCITATION’ 

3. NONLINEAR TRANSIENT DYNAMICS 
0 NEWMARX-B 

* Assumes linear load variation during each inaement 
’ With phase defined independently at each node 


NESSUS 

PROBABILISTIC ANALYSIS FEATURES 


0 THREE LEVELS OF SOPHISTICATION 

LEVa I : HOMOGENEOUS RANDOM VARIABLES 

LEVaa: stochastic variables AT THE aEMENTLEVa 
LEVa 3. STOCHASTIC INTERPOLATION WITHIN aEMENT 


0 TWO ALTERNATIVE ANALYSIS STRATEGIES; 

A. PROBABILin INTEGRATION 

B. SIMULATION (MONTE CARLO) 
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NESSUS 


EFFiciDff scamm of the 

PERTURBED STRUCTURE PROBLEM 

A. aASSlCALNONLINBARFEM PROBLEM 
o EOUILIBRIUM 

1(0) -P 

o NEWTON ITERATION METHOD 

• B^*^ ♦ vilh B^®^ ■ 0 

B. PERTITOED STRUCTURE PROWEM 
o EOUILIBRIUM 

1^(0). P 

o NEWTON-TYPE ITERATION 

with 
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NESSUS 


NUMERICAL EXAMRE 
Squire Plite. Qiinped Edges 


1 in thickness 

E - 30000 ksi 
n»03 
p * 10 psi 


10 9 10 in 
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Reducing the thickness at 1 1 nodes along one clamped edge 



Center Deflecuon (in) 

Residual on 

DiSDlacemeni 

at Iteration 


JQIEfidL 

ITERATED 




by 5t 

-0 4615 

-0.4614 

0.038 

0.001 

n/a 

by lot 

•0 4649 

-0 4644 

0.071 

0.003 

n/a 

by I5t 

-0 4669 

-04672 

0.099 

0.010 

n/a 

by 20t 

-0 4711 

-0 4710 

0.124 

0013 

0.002 

byZ3t 

-0 4740 

-0 4736 

0.144 

0.021 

0 003 


with It tolerance on the displacement residuals 
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' N85-27957 

PROBABILISTIC FIHITB ELOIBBT - VARIATIONAL THBORY* 

T. Belytsehko and Uins Kan Liu 
Horthweatarn University 
Evanston, Illinois 60201 


Traditionally, uncertainty analysis in structural mechanics has concen- 
trated on problems of an almost totally stochastic nature. Within this 
setting, even a single degree of freedom system with nonlinearities poses a 
formidable challenge and has not been solved satisfactorily. However, in 
probl»ns such as the nonlinear analysis of rocket engine components, concern 
lies more with deviations in loads from a deterministic path and in effects of 
uncertainties in material properties and boundary conditions. Thus the goal 
of this research is to provide techniques which are cost-effective and enable 
the engineer to evaluat the effect of uncertainties in complex finite elanent 
models. 

It is important to be able to treat the effects of uncertainties in a 
reasonably economical manner; standard Monte Carlo procedures are simply too 
expensive. Furthermore, the "lethods should be designed so that they can be 
incorporated into widely used finite element programs in a natural and concise 
manner. Thus, the approach should be integrable with the elemental discret- 
ization and nodal assembly procedures that characterize finite element theory 
and software implementation. 

For these reasons, embedding the probabilistic aspects in a variational 
formulatiw is a natural approach. In addition, a variational approach to 
probabilistic finite elements enables it to be incorporated within standard 
finite element methodologies. Therefore, once the procedures have been 
developed, they could easily be adapted to existing general purpose programs. 
Furthermore, the variational basis for these methods enables them to be 
adapted to a wide variety of structural elements and to provide a consistent 
basis for incorporating probabilistic features in many aspects of the 
structural problem: i.e., in the displacements, boundary conditions, body 
forces resulting from acceleration loads and any other features that cannot be 
clearly established, for example, the well known dilemma as to whether a 
shell is clamped or simply-supported at a boundary, could also be treated more 
rationally by using a probabilistic distribution for this boundary condition. 

The tasks concluded in the first year's study include the following: 

1. Theoretical Development of Probabilistic Variational Equations for 
Structural dynamics . We have developed a methodology that can embed the 
probabilistic character of the constitutive properties and leads (i.e., 
material uncertainties and load uncertainties) with a finite element 
variational approach. The corresponding probabilistic character of the 
elemental nodal forces can then be assembled into a description of the 
probabilistic distribution of the nodal forces for the complete model and the 
appropriate mean response due to these uncertainties can be determined 
efficiently. The basic concept of this approach is to incorporate the 
probabilistic distributions, as reflected in the variance, of the material 
properties and the loading to obtain the corresponding variances in the 


*Work performed under NASA grant NAG3-535. 
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elemental nodal forces of the finite element model. On the basis of the 
variance in the elemental nodal forces, the variance in the final solution 
will then be determined in the usual deterministic solution procedu'*es. 

2. The Development of EfficiOTt Numerical Algorithms for Probabilistic 
Sensitivity Displacement and Stress /^na^sis . The second step of this 
approach to probabilistic analysis is to to construct an economical method of 
obtaining the probabilistic sensitivity element rat rices which reflect the 
effects of randomness on response variables such as displacements, stresses, 
etc. The randomness is due to the preassigned probabilistic descriptions of 
the material properties and loads. While separate analyses for random 
distributions of all probabilistic response variables are an obvious 
possibility to achieving this aim, the computer cost would be formidable. It 
would require n analyses (n is the number of sample points) for each set of 
probabilistic variables. Obviously, more effective procedures must be 
devised. Relatively efficient solution procedures which have be«i developed 
here »ill be demonstrated with results for both a two-degrees of freedom 
spring-mass system and a ten-bar nonlinear syston. 

3. Integration of Methodologies into a Pilot Computer Code . The computer 
solutions obtained using this theory have been compared to the Monte Carlo 
methods and the Hermite-Gauss Quadrature integration schemes. The cost of the 
new method is substantially lower. We have demonstrated this with a ten-bar 
probabilistic nonlinear system where the random variables are the yield 
stresses. Exploitation of this characteristic in any probabilistic analysis 
r.ay perhaps offer significant savings and we will further investigate this in 
the probabilistic finite element method. 



Schematic Illustration of the Probabilistic Finite Element Method (PFEM) 
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Schematic of Probabilistic Methods 
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N85-27958 

PROBABILISTIC STRUCTURAL ANALYSIS OF SSKE TURBOPUHP BLADES - 


PROBABILISTIC GEOMETRY EFFECTS* 


Vinod K. Nagpal 
Sverdrup Technolo'jy, Inc. 
Hiddleburg Heights, Ohio 44i20 


A probabilistic study has been initiated at L«vis Research Center, "lie 
first objective of this stu<^ is tc evaluate the precisions of the geonetric 
and material properties tolerances on the structural response of turbopin^) 
blades. To complete this stu<ty» a number of important probabilistic 
variables have been identified which are conceived to affect the structural 
response of the blade. In addition, a methodology has been developed to 
statistically quantify the influence of these probabilistic variables in an 
optiniize*., way. 

The identified variables include randran geonetric and material properties 
perturbations, different loadings and a probabilistic conbination of these 
loadings. Influences of these probabilistic variables are planned to be 
quantified by evaluating the blade structural response. The structural 
response includes natural frequencies, root (maximun) stress, stage wei^L 
and tip displacements. 

Studies of the geonetric perturbations have been conducted for a flat plate 
geometry as well as for an SS^C blade oeoraetry using a special purpose code 
which uses the finite elanent approach. The geonetric perturbations which 
sinulate the natural perturbations under operathig conditions, are generated 
bi’ randomly perturbing the x, y, and z coordinates of all the nodes of the 
finite elenent mesh. Analyses of both the gecanetries indicate that the 
variances of the perturbations about given mean values have significant 
influence on the response. However, for neither of the two geometries ill 
three variances cf perturbations along x, y, and z axes showed significant 
effect. Also, the two analyses indicated that the perturbation variances 
which influence the flat plate response significantly were different from 
those of the SSME blade gecmetry. Contrary to the influence of variances, 
the means of perturbations have insignificant influence on the resporue of 
either geometry. 

Statistical mo<3els have been developed to prec3ic± the structural response 
for given perturbation means and variai.oes for both gecmecries. Since meai'.s 
arc not significant, statistical models based only o.i Ure significant 
variances have also been developed. In addition, probability distributions 
for the structural response have been developed. These cSistributions mostly 
appeared discrete. These distributions provide an assessrrent of variation 
in the structural respo*'se for a selected geometric perturbation. Also, 
these provide an estimate of the probability of getting a ce*"tain response 
for a given input. 


*Work performed under NASA Lewis task number 84-52-03. 
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Statistical tests and methods were used to <dieck the developed models. 
Ihese tests indicated that the developed models are good fits. t-Tests were 
used to identify the significant variances of perturbations. F-tests, and 
plots and autocorrelations of residuals were used to check the goodness-of- 
fit tests. 

Die preliminary conclusions frcm the studies completed so far are that the 
means of geometric perturbations are insignificant. However, the variances 
of the perturbations have significant effect. Also an an^ysis ii^cates 
that the structural response varies by up to ten percent when a test run 
with the same mean and variance of a perturbation is r^icated. 

Itie studies of randcan perturbations in material properties, different 
loadings and their probabilistic cxanbinations is underway. 


PRESENT OBJECTIVE 

CONDUCT PROBABILISTIC STRUCTURAL ANALYSIS OF 
SSflE BLADES USING PROBABILISTIC DESCRIPTION OF 

geometry 

MATERIAL PROPERTIES 
LOADING CONDITIONS 


PROBABILISTIC VARIABLES 

• :.AND0M GEOMETRIC PERTURBATIONS 

• RANDOM MATERIAL PROPERTIES PERTURBATIONS 

• LOADING 

OETERMISISTIC 

PROBABILISTIC 
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BUDE ROOT 


BLADE TIP 



NUMBER OF ELEMENTS = 80 
NUMBER OF NODES = 55 

FINITE ELEMENT MESH '^OR FLAT PLATE GEOMETRY 


SUMMARY OF STATISTICAL MODELS 
FOR SSME . BLADE. CEOMETRY 

MODEL: 

DFP. VAP. - CONSTANT + COEFF^M^ + COEFF 2 M 2 + COEFF 3 U 3 + COEFF^CT^ + C0EFF5CT2 COEFF 0 (T 3 


DEPENDENT 

VARIABLE 

SYMBOL 

CONSTANT 

COEFFICIENTS OF 

CALCULATE 

F-VALUE 

TABLE 

F-VALUE 

95Z 

Ml 

W2 

Uj 

<r. 

(T. 

■^5 

FIRST FREQ. 

'^1 

6219.9 

-2312. 

-Ikll, 

-974. 

-8716.3 


-8790. 

23.01 


SECOND FREQ. 

F^ 

c 

9923.3 

-^4517. 

-5615. 

492. 

-11094 

-6375. 

-12612. 

7.67 


THIRD FREQ. 

h 

15723.9 

-5926. 

5357. 

-2671 

-24816 

-12560. 

-36929. 

27.34 

2.26 

ROOT STRESS ! 

RS 

57039. 

19A0. 

257695. 

71492. 

11772C. 

123200. 

235795. 

2.83 


STAGE WEIGHT 

WT 

0.90731 

1 . 36A2 

1.827 

0.4323 

5,4521 

2.4369 

5.6581 

46.2 


TIP DISPL. 


-.0058 

0.0279 

0.288 

0.144 J 

BESI 

0,0572 

-.0412 

1.3C 
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SAMPLE OF 0A!A ACQUISITJGN oN 6E0M£T:L1C PEPiURdATION 
OF SSME BLADE 



GEOMETRIC PERTURBATIONS 



RESPONSE OF 

IHE BLADE 




MEANS 

blD. 

DEVIATION 

NATURAL FREQUENCIES 

ROOT 

stage 

DISPLACEMENTS 

Mk 

Ml 


Oy 

<Jz 

FIRST 

SECOND 

THIRD 

nRESS 

WEIGHT 


C.O 

C.O 0.0 

.02 

9.0 

.02 

6396.7 

8770.6 

IE882. 

^•3776. 

1.1435 

.ooo; 






6139.7 

8929.3 

16118. 

6^504. 

1.1383 

.00CS3 






5901.5 

8943. 5 

18173. 

7378F. 

1.1666 

' 0G2U7 






5952.4 

86'|R.2 

1%51. 

80025. 

1.1598 

.00815 

1 

0.0 

.C2 0.0 

.02 

0.0 

.02 

6139.1 

S3CJ.5 

163U9. 

61953. 

1.14^^7 

1 

I 

.00066 






6334.5 

9286.9 

16/71. 

64701 . 

1.1552 

. CG065 






65b‘).0 

9190... 

16229. 

70829. 

1.1451 

-.0001 






5950.9 

86C0.9 

15894. 

62124. 

1.1762 

-.UCC59 

C.O 

O.C O.I 

. 02 

0.0 

.02 

5834.9 

9184.9 

14240. 

60456. 

1.236C 

.00506 






6202.7 

8989.1: 

15530. 

56356. 

1.2338 

.00211 






5775.. 

9044.0 

15113. 

50945. 

1.2054 

.00248 






5671.1 

8878.7 

14580. 

89466. 

1.2105 

.00328 


DEFINinONS 


F - TEST 

■ ^ _ {RESIDUAL SUM OF SaUAR£S)/D.F. 
(REGRESSION SUM OF SQUARES )/D.F. 


T - T£ § I 


VALUE - KEAN 
bF. DEV. 
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PR03APILITY D’STRIBlITIOrfS FOR TEST RUN NO. ^2 


SAWPLE SIZE = 0 



FIRST ri;.T!''?AL FPFOi Errv 



' r*’f; 
STANTARr 


o-:2.c 

rfMA^ic:; :c: 


5 


00 5000 6000 7CC0 8000 9COO 


SECON? ?.VTURAL FPEOV^NC'^ 


P^o^ablllty Distributions For Test Run No. A2 


F'ean 37": 
St. rev. 33; 


Sample Size 60 





OraOfNAL PAQE iS 
OF POOR QUALITY 


ReslduoJ ulstrituUon for NaturoJ Frc,^cncy 
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U85-27959 

DYNAMIC CREEP BUCKLING: ANALYSIS OF SHELL STRUCTURES 


SUBJECTED TO TIME-DEPENDENT MECHANICAL AND THERMAL LOADING* 


G.J. Simitses, R.L. Carlson^ and R. Riff 
Georgia Institute of Technology 
Atlanta, Georgia 30332 


The objective of the present research is to develop a general mathematical 
model and solution methodc\ogies for analyzing the structural response of thin, 
metallic shell structures under large transient, cyclic, or static thermome- 
chanical loads. Among the system responses associated with th^se loads and 
conditions are thermal buckling, creep buckling, and ratcheting. Thus geo- 
metric and material nonlinearities (of high order) can be anticipated and must 
be considered in developing the mathematical model. 

A complete, true ab-initio rate theory of kinematics and kinetics for 
continu«*m and curved thin structures, without any restriction on the magnitude 
of the s rains or the deformations, was formulated. The time dependence and 
large strain behavior are incorporated through the introduction of the time 
rates of metric and curvature in two coordinate systems: fixed (spatial) and 

convected (material). The relations between the time derivative and the covar- 
iant derivative (gradient) have been developed for curved space and motion, so 
the velocity components supply the connection between the equations of motion 
and the time rates of change of the ^^'etric and curvature tensors. 

The metric tenser (time rate of change) in the convected material coordi- 
nate system is linearly decomposed into elastic and plastic parts. In this 
formulation a yield function is assumed that is dependent on the rate of change 
of stress, metric, temperature, and a set of internal variables. Moreover, a 
hypoelastic law was chosen to describe the thermoel? ^tic part of the 
deformation. 

A time- and temperature-dependent viscoplasticity model was formulated in 
this convected material system to account for finite strains and rotations. 

The history and temperature dependence were incorporated through the introduc- 
tion of internal variables. The choice of these variables, as well as their 
evolution, was motivated by phenomenological thermodynamic considerations. 

The nonisothermal elastic-viscoplastic deformation process was described 
completely by ” thermodynamic state** equations. Most investigators (in tlie area 
of viscoplasticity) employ plastic strains as state variables. Our study shows 
that, in general, use of plastic strains as state variables can lead to incon- 
sistencies with regard to thermodynamic considerations. Furthermore, the 
approach and formulation employed by all previous investigators lead to the 
condition that all the plastic work is completely dissipated. This, however, 
is in contradiction with experimental evidence, from which it emerges that part 
of the plastic work is used for producing residual stresses in the lattice. 


*Work performed under NASA grant NAG3-534. 
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which, when phenomenologically considered, causes hardening. Both limitations 
are not present in our formulation because of the inclusion of the *'thermo- 
dynamic state** equations. 
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W85-27960 

STRUCTURAL TAILORING OF SSME BLADES (VANES)* 


R. Rubinstein 
Sverdrup Technology, Inc. 
Middleburg Heights, Ohio 44130 


Hie engine blade design optimization progran STAEBL (Scructural Tailoring of 
Engine Blades) is available at the NASA Lewis canputer facility. The 
analysis capabilities of this program were extended to typical loading 
conations for SS!£ turbopui^ blades including themal and pressure loading. 
Input files for representative SSME blade designs were developed and sample 
optimization studies for these blades completed. 

The structural tailorl«ig progrem ocxnbines a general optimization package and 
^ finite element blade analysis package. The analysis package's 
capabilities include natural frequency, roaximun stress, and forced response 
ccnnputation, and fatigue life and flatter analysis. 

Optimization is performed using the "feasible directions" method. The 
current design is modified by perturbing the design variables so that the 
design constraints are satisfied white the objective function, such as blade 
weight, is reduced at the maximum rate. Hie program's geometric design 
variables include blade thickness distribution, thickness to chord ratios, 
and root chord. Special design variables are included for composite and 
hollow blades. Typical design constraints include natural frequency margins 
imposed to avoid resonances near integral multiples of critical engine 
speeds, maximum stress limits, flutter margins, and forced response 
constraints imposed to avoid fatigue failure. 

Hie ^ogram was designed for application to gas turbine fan blades. In 
order to extend the program's analysis capabilities to SSME blades, thermal 
and pressure load analyses, and the analysis of temperature dependent 
mechanical properties were added to the program. Optimization studies for 
flat plates and for compressor blades demostrated that such loads can cause 
significant increases in the wei^t of optimized designs. Pressure loads 
can greatly increase root stresses, and therefore require additional 
material at the blade root. Thermal loads influence the optimal design 
primarily by decreasing the blade's modulus of elasticity. This softening 
decreases the blade's natural frequencies and therefore requires hi^er 
blade wei^ts in order to avoid resonances with engine speed. 

SSME blade mocels provided by Lockheed and Rocket<^ne were converted to the 
it^xit format required by the structural tailoring program. Weight savings 
of approximately ten percent were achieved using optimized designs. 

Design reooirroendations cannot, of course, be based on these calculations 
alone. Input from the SSME blade designers is crucieil, especiedly to the 
formulation of constraints which correctly model the blade design 
requirements. However, the structure of the program is sufficiently 
flexible that ary further progran modifications required to accomodate SSME 
turbopuip blade design requirements will surely be feasible. 


*Work performed under NASA Lewis task number 84-52-01, 
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OBJECTIVE 


MODIfY STRUCTURTL TAILORING PROGRAM (STAEBL) FOR 
AFPliCATION ro SSME TURBOPUMP blades 


STAEBL 


OPTIMIZATI'N PROCEDURE 




INPUT 

I 

PREPROCESSOR 

I 

EVALUATE CONSTRAINTS 




/ 

GraC'^^nt of 
objective 


Not 
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Satisfied Satisfied 




\ 


Gradient of 
constraints 


ID Search 
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Ootimized Design 
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Convergence Check 
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^NOt 
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STOP 
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Feasible ^ 
design 
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STAEEL 


OWGH.'.-.L 'i: 
OF FocR Quality 


TYPICAL DESIGN VARIA^^LES TYPICAL CONSTRAINTS 


Thickness ut 5 stations 
Root Chord 

rMckness-ro-chord ratiu 


Resonance margin for order 
number and root 

(Campbell Diagram) 

Maximum roc- static srress 

T:p mode freauency margiti 

Forced response margin 
(Goodmon Diagram) 


OPTIMUM BLADE DLSiTN 
Sixth stage Compressor Blade 


INITIAL DESIGN OPTIMIZED DESIGN CONSTRAINT VALUES 


STAGE WEIGHT (LBS) 
PERCENT SPAM 
THICKNESS (IN) 
CHORD ‘IN) 
THICK/CHORD 

11. :z 

OX 50X 1002 
.29 .20 .11 
2.8 3.1 3,6 
.10 .00 .03 

7.95 

OX 5CX iOOX 
.29 A'S .08 
1.7 1.9 2.2 
.19 .0/ .09 

1 

.0" ,02 .02 

.15 .12 ,10 

CONSTRAINTS 




RESONANCE MARGINS 




MODE 1 

.23 

.05 

> .05 

MODE 2 

.21 

1.3 

> .05 

MODE 5 

.29 

1.5 

> .05 

ROOT STRESS (KSI/ 

36 

25 

< 97 
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OPTIMUM BLADE DESIGN 
Sixth Stage Compressor Blade 
EFFECT OF THERMAL AND PRESSURE LOADS 


NO THERMAL LOADS THERMAL LOADS THERMAL AND 

PRESSURE LOADS 


STAGE-WEIGHT (LBS) 
PERCENT SPAN 
THICKNESS (IN) 
CHORD (IN) 
THICK/CHCRD 

7.95 

GZ 50 % 100% 
.2A .13 .08 
1.7 1.9 2.2 
.14 .07 .04 

8.24 

0% 50% 100% 
.24 .15 .08 
1.7 1.9 2.2 
.14 ,08 .04 

9.69 

Ot 50% 100% 
.48 .08 .08 
3.3 3.6 4.2 
.14 .02 .02 

CONSTRAINTS 


1 


RESONANCE MARGINS 




MODE 1 

.05 

.06 

.05 

MODE 2 

1.27 

1.27 

1,47 

1 MODE 5 

1.65 ! 

1.C5 

1.67 

j ROOT STRESS (KSI) j 

26 j 

28 

37 


CDTIMUM BLADE DESIGN 


Lockheed SOME Blade 


INITIAL DESIGN 


O^TTMIZED DESIGN 
High root stress 


OPTIMIZED DESIGN 
Lov^ root stress 


BLADE WEIGHT (LBS) 
NUMBER OF BLADES 
PERCENT SPAN 
THICKNESS (IN) 
CHORD (IN) 
THICK/CORD 

.040 

76 

0% 50^ 100% 

.23 .15 .12 

1.0 .85 .81 

.22 .18 .14 

.026 

100 

0% 50% 100% 

.1^ .15 .06 

.79 .64 .61 

22 .23 .10 

.058 

82 

C% 50% 1C0% 

.24 .15 .11 

.96 .77 .74 

.25 .17 .15 

CONS rPAI NTS 




RESONANCE MARGINS ’ 




MODE 1 ; 


2.5 

2.9 

nODE 2 i 


4.2 

4.2 

MODE 3 


8.1 

8.6 

ROOT STRESS (KSI) 


145 

99 
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OVERVIEW OF FATIGUE t FRACTURE, AM) CONSTITUTIVE MODELING PROGRAMS AT LEWIS 


Gary R. Halford 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


Because of the absolute requirement for having safe and durable reusable 
space propulsion systems and because of the extremely high performance required 
of such systems, the means must be available for accurately calculating the 
extent of usable mission lifetime. The design of reusable space propulsion 
systems is such that many of the critical components cannot avoid having very 
short lives. These lives must be reasonably well known in advance of commit- 
ting a design to hardware because of the exceptionally high expense involved 
and the long time required for construction and certification. A minimum life 
must be attainable or the entire effort is not economically feasible. Being 
able to calculate expected lifetimes from expected loadings and structural 
analyses of the local stresses and strains enables the designer to alter 
details of a design to permit the maximizing of useful cyclic lifetime in con- 
junction with other design parameters such as performance, weight, and size. 
Knowledge of the possible design tradeoffs and their quantitative assessment 
will give rise to better designed systems. 

Hence, the objective of the NASA Lewis Research Center* s program in this 
durability area is to tf.evelop and verify cy *ic constitutive stress-strain and 
life-prediction models for materials used ir >usable space propulsion systems. 
The life-prediction models are aimed at accu ely calculating the initiation 

and early growth resistance of cyclic cracks . materials used in hot-gas-path 

components. The development of cyclic crack ropagation analysis and life- 
prediction methods by using fracture mechanics principles is being handled by 
my counterpart at the Marshall Space Flight Center, Carmelo Bianca. 

Our approach to meeting the exceptionally rigorous demands of space pro- 
pulsion systems is to start with the technology base developed in support of 
the aeronautics gas-turbine industry. Thus, life prediction and cyclic con- 
stitutive material models developed for aeronautics propulsion systems are 
being examined for deficiencies relative to reusable space propulsion systems. 

Existing models are being modified and new ones will be developed as nec- 
essary to describe the evolutionary development of fatigue damage in terms of 
phenomenological variables and micromechanistic features of the materials 
involved. Crack initiation approaches are being addressed from experimental 
and theoretical points of view. Models will be subjected to rigorous verifi- 
cation at the subcomponent level to demonstrate their applicability to the 
unique conditions of reusable space propulsion hardware. 

Vie are in the second year of the program and results are just now begin- 
ning to be realized. The early results will be presented by principal inves- 
tigators from Lewis, universities, other Government laboratories, and industry. 


rRECKDlNO r\C,V BLANK NOT FU.MED 
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W85-27961 

INTERACTION OF HIGH-CYCLE AND LOW-CYCLE FATIGUE OF HAYNES 188 ALLOY AT 1400 


Peter T. Bizon, Dan J. Thoma» and Gary R. Halford 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


The interaction of low-cycle fatigue (LCF) and high-cycle fatigue (HCF) 
was evaluated at the NASA Lewis Research Center on Haynes 188 alloy at 1400 ""F. 
Completely reversed, axial-load, strain-controlled fatigue tests were performed 
to determine the baseline data for this study. Additional specimens for inter- 
action tests were cycled first at a high strain range for various small por- 
tions of expected LCF life followed by a step change to a low strain range to 
failure in HCF. Failure was defined as complete specimen separation. The 
resultant lives varied between 10 and 5000 cycles for the low-cycle fatigue 
tests and between 4500 and 3 million for the high-cycle fatigue tests. For the 
interaction tests the low-cycle-life portion ranged from 30 and 1000 applied 
cycles while the high-frequency life ranged from 300 and 300 000 cycles to 
failure. The step change results showed a significant nonlinear interaction 
in expected life. Application of a small part of the LCF life drastically 
decreased the available HCF life as compared with what would have been expected 
by the classical linear damage rule (LDR) . 

The nonlinear cumulative damage rule proposed by Hanson and Halford in 
1981 and referred to as the ’’damage curve approach (DCA)” predicted the trends 
of the results. However, the presently observed interactions were more severe 
than anticipated, implying that the single universalized constant in the DCA 
has a different value than found for previous materials tested. A fracto- 
graphic study of failed specimens showed that the fracture surfaces from inter- 
action tests appeared to be essentially the same as surfaces from specimens 
tested in HCF only: the LCF portion of the cycling was not distinguishable on 
the interaction fracture surfaces. 
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NOMINAL PROPERTIES OF HAYNES 188 ALLOY 
CO M P OS IT ION (wt %) 


COBAU 

39^ 

IRON 

3% MAX. 

NICKEL 

22% 

MANGANESE 

1.25% MAX. 

CHROMIUM 

22% 

SILICON 

0.4% MAX. 

TUNGSTEN 

14% 



CARBON 

0.1% 



LANTHANUM 

0.08% 




A T 1400 °F (TEST TEMPERATURE! 

YIELD SIRENGIH (a 2% OFFSETI 42000 psi 

ULTIMATE TLNSIIE STRENGTH 92000 psi 

TENSILE ELONGATION 43 % 

V-2490 


HOURGLASS SPECIMEN GEOMETRY 

DIMENSIONS IN INCHES 



V-2485 

CS-84-450c 
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TEST PARAMETERS 


COMPLETELY REVERSED, AXIAL LOAD, STRAIM-CONTROl., FATIGUE TEST 

TEMPERATURE 1400 

HEATING METHOD SELF RESISTANCE 

FREQUENCIES a4 AND 29 Hz 

LOAD TRAIN ALIGNMENT “WlOODS-METAL SYSTEM" 

MEASUREMENT: 

LOAD FOUR-ARM STRAIN GAGE LOAD CELL 

STRAIN DIAMETRAL EXTENSOMETER USING LVDT 

TEMPERATURE OPTICAL DISAPPEARING-V\IIRE TYPE 

PYROMETER 


V-2489 


CS-84*4302 


PROGRAM SCOPE 


FATIGUE TESTS: 

NUMBER OF TESTS 

FREQUENCY (Hz) 

LIFE (CYCLES) 

9 

a4 

13 -4918 

n 

29 

4500 —3262300 

10 

a 4— 29 

33-925^300-310800 

METALLOGRAPHIC STUDY: 

TYPE OF TEST 


LIFE (CYCLES) 

LCF ONLY 


2142 

HCF ONLY 


3262 300 

COMBINED LCFIHCF 


51—76000 


V-2488 


CS-ft ,^4903 
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SSME STRUCTURAL DURABILITY 


DAMAGE CURVE APPROACH 

TWO LEVELS OF LOADING 
A N2 


fNl 

1- \Hl] Ha 

Hi 


-INITIAL LIFE FRACTION 


REMAINING LIFE FRACTION 
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OftiCNNAL PAGE tS 
OF POOR QUALITY 


SSME STRUCTURAL DURABILITY 

DAMAGE CURVE APPROACH(DCA) 

TWO LIFE LEVELS, N,^AND 



(a) Ni = N,^ . N2= Nhi (b) Ni= N^; . N2=N,„ 


BASELINE FATIGUE DATA 

HAYNES 188 ALLOY AT 1400 ^ 


TEST 

FREQUENCY, 

Hz 
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RBEXAMINATIOa OF CUMULATIVE FATIGUE DAMAGE LAWS 
Gary R. Halford 
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S.S. Manson 
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Cleveland, Ohio 44106 


Treatment of accumulated fatigue damage in materials and structures sub- 
jected to a history of nonsinq>le repetitive loadings has received a large 
amount of attention in recent years. The subject has been a popular one for 
more than half a century since Palmgren (ref. 1) first suggested what is now 
known as a "linear damage rule." After the same rule was later independently 
proposed by Langer (ref. 2} and Miner (ref. 3), it vas recognized that, 
although the method has much merit for the simple treatment of complex loading 
history, it can result in unconservative predictions of material and structural 
behavior. An intense flurry of activity followed in the pursuit of alternative 
methods of analysis that would predict behavior more accurately. So many meth- 
ods were introduced that it became necessary periodically to prepare review 
papers placing all the new methods into perspective, among them the ones by 
Newmark (ref. 4), Kaechele (ref. 5), Hanson (ref. 6), O’Neill (ref. 7), Schive 
(ref. 8), and Laflen and Cook (ref. 9). Schive's study lists nearly 200 ref- 
erences while reviewing approaches to treat the problem and the numerous exper- 
imental programs that are relevant to the evaluation of the concepts. No truly 
con^rehensive review has appeared of the considerable effort that has been made 
in the pest dozen years since Schive’s publication. 

The current presentation is based on a few of the results contained in an 
extensive paper being prepared on the subject by Manson and Halford (ref. 10). 
Our primary goal is to review our own direct experience with this subject, to 
briefly mention some recent progress, and to express our current integrated 
view regarding the state of the art as it applies to our own efforts. Only 
incidentally shall we review the more recently proposed cumulative damage life 
prediction methods. 

The double linear damage rule (DLDR), which has evolved in our laborato- 
ries over the past 20 years (refs. 11 to 13), is reexamined with the intent of 
improving its accuracy and applicability to engineering problems. Modifica- 
tions are introduced to the analytical formulation to achieve greater compati- 
bility between the DLDR and the so-called damage curve approach (ref. 13), 
which is an alternative continuous representation of the DLDR. For our present 
purposes we shall demonstrate the new formulation by using the recent HCF/LCF 
results generated by Bizon et al. (ref. 14) and presen’ ed at this conference. 
The constant in the new formulation are derived directly from the new data on 
the cobalt-base alloy Haynes 188. Once calibrated, the equations are used to 
predict lifetimes under complex, yet more realistic, simulated service 
conditions. 
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DOUBLE D^. iGE CURVE APPROACH(DDCA) 
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^ — Small Integer as Required 


Previous Work 

oC = o.as ; yS = o.40 

Haynes 188 at 1400 F 

oC = 0-35 ; 3 = O.bO 
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CORRELATION OF HCF/LCF INTERACTION BY THE 
DOUBLE DAMAGE CURVE APPROACH(DDCA) 



DOUBLE LINEAR DAMAGE RULE(DLDR) 

• 0 r- 
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PHASE I ("Initiation”) 



PHASE II ("Propagation”) 



SIMULATED LOADIJ^G CASE(SSME LOX Post) 

• HCF Freq. = lOuunz 

• Mission Duration = SOOsec. 

• LCF Life = 500 Cycles 

• HCF Life = 50,000,000 Cycles 

Mission Consists of 1 LCF Cycle at N^= 500 
and 500,000 HCF Cycles at 50,000,000 

Based on Miner Linear Damage Rule, 
n/N,=: 0.005 and n/Nj= 0.01 
0.005-!- 0.01= 0.015 
No. Missions = 1/0.015 - 67 
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Based on DLDR for Haynes 188 at 1400 F 
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Fuel turbopurop components for reusable space propulsion systems experience 
extreme gas temperature and pressure environments. These operating conditions 
subject the high-pressure-stage turbine nozzles and blades to severe thermal 
transients that result in large inelastic strains and rapid crack initiation. 

To attain even the short lives required for these systems, anisotropic turbine 
blading alloys must be used. Assessing or improving the durability of these 
components is contingent on accurate knowledge of the stress-strain history at 
the critical location for crack initiation. 

Nonlinear finite-element analysis techniques have become available in 
recent years for calculating inelastic structural response under cyclic load- 
ing. These methods are based on classical incremental plasticity theory with 
uncoupled creep constitutive models. Many of the nonlinear finite-element 
computer codes such as MARC (ref. 1) have the capability of handling materials 
with anisotropic properties. However, these codes are usually too costly and 
time consuming to use in the early design stages for aerospace applications. 
Costs are further increased by the geometrical complexity of high-pressure 
turbine blades, which necessitates three-dimensional analyses and sometimes 
substructuring for accurate solutions. Simplified procedures for more econom- 
ically representing structural response under cyclic loading have been devel- 
oped (refs. 2 to 4) to improve the design of engine hot-path components such 
as turbine blades. 

This paper addresses the problems of calculating the structural response 
of high-tamperature space propulsion components such as turbine blades for the 
fuel turbopump. The first high-pressure-stage fuel turbine blade (HPFTB) in 
the liquid-hydrogen turbopump of the space shuttle main engine (SSME) was 
selected for this study. In the past these blades have cracked in the blade 
shank region and at the airfoil leading edge adjacent to the platform. To 
achieve the necessary durability, these blades are currently being cast by 
using directional solidification. Single-crystal alloys are also being inves- 
tigated for future SSME applications. The study evaluated the utility of 
advanced structural analysis methods in assessing the low-cycle fatigue lives 
of these anisotropic components. 

The turbine blade airfoil of the high-pressure stage of the SSME fuel 
turbopump was analyzed because it has a history of rapid crack initiation. 

These uncooled airfoils have a span length of 0.33 cm and a span-to-chord width 
aspect ratio of approximately unity. The blades are cast from directionally 
solidified MAR-M 246 alloy. Temperature-dependent properties for this alloy 
were mainly provided by the Rocketdyne Division of Rockwell International Cor- 
poration. Material elastic properties are summarized in table I. Mean thermal 
coefficient of expansion data were converted to instantaneous values for MARC 
input. Poisson’s ratio values of 0.143 and 0.391 were used for the longitudi- 
nal and transverse directions, respectively, and were a<^sumed to be constant 
with temperature. Longitudinal stress-strain properties, summarized in 
table II, were used for the elastic-plastic region. A single- crystal MAR-M 
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246 alloy is also being considered for turbine blades in future SSME 
applications. 

Cracking has occurred during service at the airfoil base near the leading 
edge and in the blade root shank area. These cracks apparently initiate during 
the first few mission cycles because of the severe thermal transients and are 
propagated by vibratory excitation. Since the primary purpose of this study 
was to compare nonlinear finite-element and simplified analytical methods » the 
blade root and platform were excluded from the analysis to reduce the size of 
the problem and therefore the computing time. 

The mission used for the analysis is shown in figure 1 in terms of turbine 
inlet temperature, gas pressure, and blade rotational speed. This cycle is 
applicable to a factory test of the engine; it is also reasonably representa- 
tive of a flight mission except for the foreshortened steady-state operating 
time. The major factors inducing rapid cracking are the transient thermal 
stresses and the inelastic strains caused by the sharp acceleration and 
deceleration transients. 

Rudimentary transient and steady-state three-dimensional heat transfer 
analyses have been conducted using the MARC code. Film coefficients were 
obtained from preliminary information supplied by Rocketdyne. The gas sink 
temperature was assumed to be constant around the airfoil suL.ace for each time 
step. Colder boundary conditions were assumed at the airfoil base to simulate 
the effects of the cooling of the blade-to-disk attachment region by the 
liquid-hydrogen fuel. Calculated meta? temperatures at the leading edge at 
midspan and at the crack initiation site at the base of the airfoil (critical 
location) are presented in figure 2 as a function of elapsed time during the 
cycle. The assumed gas sink temperature around the airfoil is also p^^otted. 

Of particular note is that the leading-edge temperature is lower at the airfoil 
base than at midspan throughout the cycle. This seems reasonable because of 
the cooling of the blade-to-disk attachment region by the liquid-hydrogen fuel. 
The lower airfoil base temperatures induce tensile thermal stresses at the 
critical leading-edge location that are additive to the centrifugal stresses. 

Elastic-plastic analyses are being conducted for the HPFTB airfoil with 
both the MARC code and a simplified analytical procedure .developed at the NASA 
Lewis Research Center. The MARC analysis used the anisotropic material prop- 
erties, while the simplified analy^iis used only Young’s modulus and Poisson’s 
ratio for the longitudinal (spanwise) direction. 

A three-dimensional finite-element model of the airfoil (fig. 3) was con- 
structed of eight-node isoparametric elements. The model consisted of 360 
elements with 576 nodes and 1661 unsuppressed degrees of freedom. This model 
was a shortened version ot a finite-element model created by Lockheed (ref. 5) 
for a NASTRAN elastic analysis of the HPFTB blade at rated power level; this 
work was done under contract to the NASA Marshall Space Flight Center. The 
main difference between the NASA Lewis and Lockheed models was that the blade 
base and most of the platform were omitted for the MARC nonlineat analysis to 
reduce the computing time and to run the problem in-core on the Cray computer 
system at Lewis. Boundary conditions were applied to constrain all nodes at 
the base of the model to lie on a plane. Additional boundary conditions were 
imposed to prevent rigid-body motion in this plane. 
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The MARC code has been used extensively at Lewis for inelastic analyses 
of aircraft turbine blades and combustor liners and of space power components. 
In conducting a cyclic analysis the loading history is divided into a series 
of incremental load steps » which are sequentially analyzed. The elasticity 
algorithm is based on a tangent stiffness approach in which the stiffness 
matrix is reformulated and reassembled for every plastic 1^ increment. The 
incremental loads are modified by residual load correction « actors to ensure 
that the solution does not drift from a state of equilibrium. Convergence for 
the iterative plasticity analysis is indicated when the strain ener^^y used in 
assembling the stiffr^ess matrix approximately equals the energy result- 

ing from the incremental solution. 

The mission cycle was subdivided into 129 increments. Incremental loading 
included centrifugal and gas pressure loads and metal temperature distributions 
as calculated from the heat transfer analysis. Approximately 1 milii^n words 
of core storage on the Cray system were needed to run the problem. Each cycle 
of analysis required about 15 hours of computing time. 

The directionality of the elastic material properties causes anisotropic 
constraints. Lekhnitskii (ref. 6) has derived the genera.Hzed elastic strain 
equations for an anisotropic body with a transver^ie plane of isotropy Matrix 
inversion of these equations to solve for the stresses results in the 
relationship 
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= E/E* and a 

= nE'/(l + v)(l - V 

- 2nv*2) 

Here 

E* 

G* , and 

V* 


denote Young *s modulus » the shear modulus, and Poisson* s ratio, respectively, 
for the longitudinal or span direction while E, G, and v denote these 
constants with respect to any direction in the transverse plane of isotropy. 
This anisotropic stress-strain law was incorporated i** the MARC user 
subroutine HOOKLW. 

Plastic strain calculations were based on incremental plasticity theory 
using the von Mises yield criterion, the normality flow rule, and a kinematic 
hardening model. The material elastic-plastic behavior was specified by the 
yield strengths and work-hardening properties in the longitudinal direction; 
transverse properties were not available. Creep analyses are not being per- 
formed at the present time because of inadequate kno>'ledge of the creep char- 
acteristics for the anisotropic blade material. 

At the beginning of this study a steady-state MARC elastic-plastic anal- 
ysis of the airfoil was ccaducted for rated power level (RPL) . The temperature 
input was obtained from a NASTRAN steady-state heat transfer analysis performed 
by Lockheed. In contrast to the large temperature gradients calculated from 
the NASA heat transf r analysis, Lockheed had a maximum temperature variation 
in the airfoil of oni/ 35 deg C. Effective stret*:s distributions, as calculated 
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dt the Gaussian integration points closest to the suction and pressure sur- 
faces, are presented in figure 4. Since the airfoil was essentially at an 
isothermal state, the stresses primarily reflect the centrifugal and gas pres- 
sure loading and were entirely in the elastic range The maximr"' effective 
stress was only 40 percent of the yield strength. It is apparetiu that ^he 
analysis would not show a low-cycle fatigue problem if transient thermal 
effects were not considered. A MARC nonlinear finite- element analysis for the 
complete mission cycle was begun but has not been completed. 

This work is still in a preliminary stage. A transient heat transfer 
analysis of the HPFTB airfoil has been completed. Tiis analysis is rudimentary 
because of the inadequate definition of the heat transfer coef f icis^nts and 
omission of the cooled shank region in the finite- element model. However, it 
is sufficient for exercising and evaluating the structural analysis methods. 

A MARC cyclic stress-strain analysis has beer started. 

The simplified procedure will probably require further development for 
application to SSME turbopump blade problems - not only because of tne use of 
anisotropic materials, but also because of the nonproportionalit> of the 
thermal and mecnanical loading during the cycle. Experience has indicated that 
effective stresses and strains are not suitable criteria for relating multi- 
axial, nonproportiona. stress-strain states to uniaxial material properties. 
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TAdLE 1. - PHYSICAL PROPERTIES OF DIRECflOtiALLY 
SOLIDIFIED MAR-M 246 


iTenqiercture , 
*C 

1 

Modulus of elasticity, 

j 

Kean coefficient of 
thermal expansion, 
%/*C 

Longitudinal 

Transverse 

21 

131 

183 


93 

128 

179 


204 

125 

173 


316 

124 

173 

.00133 

427 

119 

166 


538 

114 

162 


649 1 

109 1 

156 


760 

103 

149 


871 

97 

142 

.00160 


TABLE II. - STRAIN PROPERTIES 
(LONGITUDINAL) OF DIRECTIONALLY 
SOLIDIFIED HAR-M 246 


Plastic strain, 
percent 

Stress, MPa, 

at - 

21 “C 

649 “C 


0.1 

800 

800 

875 

.2 

830 

855 

930 

.4 

850 

895 

965 

.6 

855 

930 


.8 i 

865 

945 

975 

1.0 j 

870 


980 
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Figure 2. - Airfoil temperature cycle. 



Figure 3. • Airfoil finite element 
model. 
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STRESS. 

MPa 

1 39.4 

2 77.9 
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8 184.1 
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(b) Pressure side. 


figure - Airfoil effective stress distributior at RPL. 
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PUA 1480 is a potential candidate material for use in the high-pressure 
fuel turbine blade of the space shuttle main engine. As an engine material it 
will be subjected to high-cycle fatigue loading superimposed on a high mean 
stress due to combined centrifugal and thermal loadings. The present paper 
describes the results obtained in an ongoing program at the Argonne National 
Laboratory, sponsored by NASA Lewis, to determine the effects of a high mean 
stress on the high-cycle fatigue behavior of this material. 

Straight-gauge high-cycle fatigue specimens, 0.2 inch in diameter and with 
the specimen axis in the [001] direction, were supplied by NASA Lewis. The 
nominal room temperature yield and ultimate strength of the material were 146 
and 154 ksi, respectively. Each specimen was polished with 1-ym diamond 
paste prior to testing. However, the surface of each specimen contained many 
pores, some of which were as large as 50 yro. Testing was carried out at room 
temperature in the laboratory air nvironment. Future tests will be conducted 
at an elevated temperature in an inert environment. 

In the initial tests, specimens were subjected to axial-strain-controlled 
cycles. However, very little cyclic plasticity was observed. For example, 
figure 1 shows the hysteresis loop at the tenth cycle of a 0 to 0.8 percent 
strain cycling test. The hysteresis loop has almost zero width, and very 
little cyclic mean stress relaxation occurred during the test. Since the 
majority of the tests were conducted at much lower strain ranges, where the 
hysteresis loops have no measurable width, it was decided to switch to 
load-controlled testing at a frequency of 20 Hz. 

Figure 2 summarizes all of the tests run to date. Tests with mean 
stresses of 0, 60, and 120 ksi were carried out. The zero mean stress test 
data fall on a linear extrapolation of low-cycle fatigue test data for the same 
material at 1600 ®F in a high-pressure hydrogen environment (ref. 1). Although 
the data are limited, the material seems to exhibit an endurance limit that is 
weakly dependent on the mean stress. Figure 3, for example, shows a Goodman 
diagram for a lif 2 of 10^ cycles. It is evident that the high-cycle fatigue 
strength cf the material is quite resistant to mean-stress effects. 

Several specimens were analyzed by scanning electron and optical micro- 
scopy after fracture. Each specimen contained a large number of micropores 
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(fig. 4, top), and crack initiation always occurred from one of these micro- 
pores. Figure 5 shows crack initiation from a micropore in a specimen (#15-3) 
that was cycled at a stress range of 120 ksi with zero mean stress. A similar 
mode of crack initiation was also observed in a specimen (#112-3) that was 
cycled at a stress range of 38 ksi with a mean stress of 120 ksi (fig. 6). 

Crack propagation was crystallographic during most of the life; it apparently 
occurred along the {111} type of plane. In both specimens #112-3 and 
#15-3, the crack that led to failure was initiated b/ a crack on a different 
slip plane (figs. 4 and 5). Crack growth also appeared to preferentially fol- 
low microporosities (fig. 7). (The distinct slip band marks along the {111} 
type of plane in the figure are for a test with a rather large (0.8 percent) 
strain range. Such slip bands were not observed in the low-strain- range, 
high-cycle fatigue test specimens.) 

Several slight differences were observed between specimen with zero and 
high mean stress. The specimen with high mean stress (#112-3) appeared to have 
cracks initiated at several sites, whereas the zero-mean-stress specimen 
(#15-3, fig. 5) had a single crack origin. The high-mean-stress specimen also 
had a rougher, more textured fracture surface, which may be indicative of a 
higher crack propagation rate. It also contained several *’steps** on the frac- 
ture surface (fig. 4, bottom). These steps were probably caused either by the 
linking of two parallel cracks on the same type of octahedral plane or by the 
intersection of cracks on nonparallel octahedral planes. 
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gure 1. - Hysteresis loop at 10th cycle for PWA 1489 specimen subjected to cycling between 0 and 
0.8 percent strain range at room temperature. 
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Figure 3. - Goodnwn diagram for PWA 1480 for lives of 6000 and 10^ cycles at room temperature. 
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figure 4. ~ Optical micrograph (top) illustrates typical micropore density in speciinen #15-3 

(unetcheoj . Scanning electron micrograph (bottom) shows a "stepped** fracture surface in specimen 
f 112-3, which was cycled at stress range of 33 ksi with mean stress of 120 ksi at room temperature. 
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Figure 5. - Fracture morphology in specimen #15-3, which was cycled at stress range of 120 ksi with 
zero mean stress at room temperature* The micropore initiated two cracks on different planes; on^ 
of these cracks led to faUure of the specimen. 
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CONSTITUTIVB CYCLIC DEFORMATION BEHAVIOR IN SINGLE-CRYSTAL AND 
DIRECTIONALLY SOLIDIFIED SSMT HIGH-PRESSURE FUEL 


1URBOPUMP AIRFOIL MATERIALS 


Walter W. Milligan, Eric S. Huron, and Stephen D. Antolovich 
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The major goal of the project is to correlate mechanical properties with 
microstructural deformation behavior and to develop models for constitutive 
response under a variety of monotonic and cyclic loading cycles, teroperaturer , 
strain levels, strain rates, and environments. The main areas of research to 
date have been the development of techniques for sample preparation and micro- 
structural characterization. In addition, mechanical tr^ting is being carried 
out and results will be reported. 

Two alloys are being studied as candidate SSHE turbine blade materials. 

"lie first is PWA 1480, which is a single crystal alloy whose nominal composi- 
tion is reported in table I. As documented in figure 1, it is a two-phase 
y/y* alloy, with the size averaging 0.75 ym after heat treatment. 

The dislocation dw^^lty in the heat-treated material is very low. The alloy 
contains about 2 vol Z, microporosity (see abstract by S. Majumd-^r et al.). 

The alloy also co:i interdendritic y/y* eutectic. Several deforma- 

tion microstructures are sho\T. in figure 2. 

The second alloy being studied is D.S. Mar-M 246 + Hf. This is a direc- 
tionally solidified material, and its nominal composition is also reported in 
table I. Lihe PWA 1480, it is a two-phase y/y* alloy, with the y’ 
size averaging 0.5 ym (fig. 3). The precipitates contain a network of inter- 
facial dislocations (fig. 3), an indication of a large mismatch with 
matrix It also contains y/y* eutectic areas at the grain boundaries 
(fig. 4), large carbides (probably of the MC type) within the grains (fig. 4), 
fine carbides (probably M^C or M 23 C^) in the grain boundaries (fig. 4), 
and submicroscopic carbides in the eutectic areas (fig. 3). This material 
contains very little microporositj' . 

The major areas of interest for the two materials will be ^ilightly dif- 
ferent. The single-crystal alloy lends itself well to fundamental deformation 
studies, since resolved shea^* stresses on slip planes are all known and only 
one grain is present. On the other hand, tVe D.S. material presents an excel- 
lent opportunity to study the effects of I'lightly misaligned grains on defor- 
mation behavior. 

The two materials will be studied by using approximately the same test 
matrix, so a good degree of direct comparison will also be possible. In addl- 
t*c. to contributing to fundamental understanding, the results of this study 
ould have profound implications for design inasmuch as constitutive models 
(based on experimental observations) will be developed. 
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TABLE I. - ALLOY COMPOSITIONS 


Element 

PWA 1480 

D.S. 


1 

Har-M 246 


1 

Composition, %it % | 

Cr 

10 

9 

Co 

5 

10 

A1 

5 

6 

Ti 

1.3 

1.5 

W 

4 

10 

Mo 


2.5 

Ta 

12 

1.5 

Hf 


1 

C 


0.15 

B 


0.015 

Zr 


0.05 

Ni 

Bal. 

Bal. 
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iEM micrographs of deformation in PWA 1480. 

(a) Slip band in the as-heat treated condition 
<111> reflection, 70.000X. (b) Extremely high 

dislocation density after tensile failure at 
UOO^F. The majority of dislocations appear to 
e in the matrix. <111> reflection, 27,OOOX. 
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TEM micrographs of as-heat treated DS Mar-M 246 + Hf. 
(a) Bright field micrograph showing intcrfacial 
dislocation networks and slacking faults. <200> reflec 
20,000X. (b) Dark field of the same area, slightly 

translated.* <300> y' superlattice reflection, 20,00OX 
(c) Bright field in a y/y' eutectic area. <200> 
reflection, lO.OOOX. (d) DarK field of the same area 
as ”c", showing coherent carbides in the matrix. 
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N85-27966 

H16H-THERMAL-GRAD1ENT SUPBRALLOY CRYSTAL GROWTH 

David D. Pearson. Donald L. Anton, and Anthony F. Giamei 
United Technologies Research Center 
Bast Hartford Connecticut, 06108 


Single, tool] -oriented crystals of PWA 1480 were processed in alumina/ 
silica shell molds in a laboratory high-gradient furnace. The furnace employs 
a graphite-resistance-heated element, a radiation baffle, and a water-cooled 
radiation trap below the baffle. All crystals were grown in vacuum (10 torr) 
and all heat transfer was radiative. The element is constructed with a vari- 
able cross section that is tapered just above the baffle to maximize heat input 
and therefore thermal gradient. A maximum alloy temperature of 1600 "^C was 
used. A thermal gradient of 130 deg C/cm was recorded at 1370 ®C just above 
the solidus of the PWA 1480 alloy. Crystal bars with 14.4- and 17.5-mm diam- 
eters were grown in alumina/silica shell molds similar to those used for pro- 
ducing gas turbine hardware. Each crystal was started from a 1.6-mm pencil 
seed at a rate of 76 mm/hr and slowly accelerated to a rate of 200 mm/hr under 
computer control. Volume percent porosity ana average pore size were measured 
as functions of distance in representative bars. Low-cycle fatigue (LCF) 
behavior was determined at UTRC, and stress-rupture properties were determined 
at the NASA Lewis Research Center. UTRC will discuss the crystal growth, 
porosity measurements, and LCF behavior. 
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N85 -27967 

ROLE OF HYDROGEN AND STRUCTURE ON CMSX-2 

C.L. Baker, J. Chene, W. Kroiop, H. Pinczolitis, 

S.H. Bernstein, and J.C. Williams 
Camegie-Mellon University 
Pittsburgh, Pennsylvania 15213 


Hydrogen environments have previously been shown to be deleterious to the 
properties of nickel-base superalloy single crystals. Tne goal of this study 
is to characterize and understand the detailed effect of hydrogen on the ten- 
sile and fatigue behavior of single crystals of CMSX-2. 

Specimens were solut ionized for 3 h at 1315 air cooled, and then aged 

in one of two ways: 980 ®C/5 h/air cool + 850 ®C/48 h/air cool (standard heat 
treatment), or 1050 ®C/36 h/air cool ^ 850 ®C/ 48 h/air cool (ONERA heat 
treatment) . 

The resultant macro- and microstructure included a pronounced dendritic 
structure with numerous low angle boundaries in the interdendritic region and 
on the average a high level of porosity with a large variation in size and 
distribution. The roicrostructure consisted of 0.5-ym cuboidal y’ pre- 
cipitates separated by thin continuous regions of matrix. The ONERA heat 
treatment led to slightly larger, more aligned, and more cuboidal y’ pre- 
cipitates - a structure known to improve the creep strength of this alloy. 

Hydrogen permeability and diffusivity in this material was found to be 
very low at room temperature. To enhance the extent of hy'^rogen penetration 
and embrittlement, hydrogen was introduced by cathodic charging at 150 ®C in 
molten salt. This was expected to increase the permeability by 2 or 3 orders 
of magnitude. This treatment produced a 50-pm-deep embrittled layer in 
charged tensile and fatigue specimens thac explicitly affected the fracture 
behavior. Knoop roicrohardness measurements of the hydrogenated surface showed 
a relative increase of about 25 percent, attributed to the lattice distortions 
associated with the steep hydrogen concentration gradient near the surface and 
consistent with the fracture appearance. 

Hydrogen concentration studies on thin samples yielded values signifi- 
cantly higher than found for similar charging conditions in nickel and stain- 
less steels. This implied a high effective hydrogen fugacity in molten salt 
and a strong trapping effect primarily at voids and y/y* interfaces. The 
existence of strong trapping at voids was inferred by the continued existence 
of an embrittlement effect and a large residual hydrogen concentration after 
desorption at 150 ®C. Trapping at y/y* interfaces as suggested by the 
appearance of these interfaces on the fract tre surfa^.:^. 

The tensile mechanical properties were found to depend on hydrogen con- 
centration. Hydrogenation in molten salts enhr.ncc<^ hydrogen* s effect on 
mechanical properties, particularly a significant decrease in strength and 
plasticity, with the greatest effect being decreased elongation, presumably an 
effect of hydrogen on crack initiation. Scatte'*' in the tensile ductility of 
both charged and uncharged smooth tensile samples was observed, probably 
because of voids and slight differences in crystal orientation. No recovery 
of the mechanical properties was observed after aging for 7 days at room tem- 
perature, in agreement with t >e ver^ ow kinetics of hydrogen desorption at 
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20 ^C. After desorption at 150 the remaining embrittlement was associated 
with hydrogen trapping in voids. Unlike previous studies, no influence of heat 
treatment on the extent of hydrogen embrittlement was seen in smooth tensile 
samples for thet.a charging and testing conditions. 

SEM observations of the brittle subsurface zone revealed a strong corre- 
lation between crack initiation sites and voids just beneath the surface. 
Preferential decohesion around the voids may have resulted from a hydrogen 
pressure buildup in such voids, while the very flat rupture features of the 
brittle zone indicated that the crack propagation was perpendicular to the 
stress axis, suggesting a preferential decohesion along the y/y* inter- 
faces. In contrast uncharged samples showed a ductile tearing of the matrix 
throughout the cross section. The presence of numerous superficial transverse 
cracks in the brittle layer suggests that while cracks can initiate very easily 
in the outer zone, their propagation in the bulk is more difficult. Thus, a 
more dramatic loss ^n tensile ductility was not observed, particularly for 
smooth specimens, since the amount of "’"ogenated material was small with 
respect to the total volxime of the sample. 

Fatigue tests at constant plastic strain amplitude were performed on sam- 
ples with two different geometries: smooth solid specimens and hollow samples 

with similar external shapes. Similar to uniaxial testing, the accumulated 
strain to fracture in the hydrogen-charged samples was lower than in the 
uncharged ones; however, the ratio of the total strain (charged versus 
uncharged) was considerably lower in fatigue than in uniaxial tension. The 
decrease in the UTS of the charged sample is similar to uniaxial testing: 
about 10 percent lower than the uncharged ones. 

The hollow and solid samples were compared to better understand the 
influence of increased charged volume on the fatigue behavior. As expected, 
the fatigue lifetime of the hollow samples was significantly shorter than that 
of the solid ones. In fact, the numbers of cycles to failure was about two 
times smaller while the ratio of the hydrogen-affected volume in the hollow 
versus solid specimens was approximate) y 4, suggesting a scaling relationship 
between the two. The fracture appearance was also similar to those of the 
tensile samples, especially the enhanced cracking around voids. 

CONCLUSIONS 

1. Sig’ ^ficant hydrogen embrittlement of CMSX-2 can occur in high~fugacity 
cathodic environments. 

2. Hydrogen accelerates the ease of crack initiation in both tension and 
fatigue. 

3. Fatigue testing magnifies the extent of hydrogen embrittlement compared 
to tension. 

4. Hydrogen embrittlement is enhanced by the presence of internal voids, 
which act as trap<? for hydrogen. 

5. No effect of heat treatment on catbodicaliy included hydrogen embrit- 
tlement of smooth tensile samples has been observed. 
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6. The dendritic structure deflects the crack path in hydrogenated 
specimens . 

7. High concentrations of hydrogen affect the work-hardening behavior of 
CMSX-2 . 

8. Permeability of hydrogen is lo^ in CHSX-2. 
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OVERVIEW OF INSTRUMENTATION PROGRAM 


William C. Nieberding 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


A great deal of emphasis in the SSME durability program is being placed 
on the development of analytic codes for designing components and predicting 
their expected life* Before any faith can be placed in a code, there must be 
a series of benchmark experiments specifically designed and instrumeated to 
verify the operating conditions that have been used as input to the code and 
to determine whether the code does indeed predict the correct results. 

The instrumentation portion of the durability program is specifically 
aimed at developing those sensors and instrument systems necessary to verify 
the codes. In most cases this calls for developing sensors and systems that 
are minimally intrusive so as not to perturb the conditions being measured. 

In addition, since the codes tend to work with high- resolution maps of 5 , :h 
parameters as flow, strain, and temperature, the instrument systems must pro- 
duce data that are very high in both spatial and temporal resolution. 

Much effort has already been put into developing such instrumentation for 
aircraft engines. This program is designed to extend this work to include the 
operating conditions and materials of the SSME. Indeed, some of these condi- 
tions arc far more stringent than in aircraft engines, but other conditions 
are less stringent. 

In this session you will hear four papers on the development of instrument 
sensors and systems aimed primarily at verifying the analytic codes being 
developed in the other parts of the durability program. 


181 




N85-27968 

THIH-PILH SENSORS FOR SPACE PROPULSION TECHNOLOGY 


Walter S. Kim and David R. England 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


Specific needs in the space shuttle main engine technology are measure- 
ments of temperature, heat flux, and strain on the components to assess their 
structural integrity and durability. SSME components such as the turbine 
blades of the high-pressure fuel turbopump are subjected to rapid and extreme 
thermal transients that contribute to blade cracking and rubsequent failure 
(ref, 1 ). 

The objective of this work is to develop thin-film sensors for SSME com- 
ponents. Thin-film sensors have been developed for aircraft gas turbine 
engines and are in use for temperature measurement on turbine blades to 1800 ""F 
(fig. 1). The technology established for aircraft gas turbine engines will be 
adopted to the materials and environment encountered in the SSME. Specific 
goals are to expand the existing in-house thin-film sensor technology, to con- 
tinue developing improved sensor processing techniques via grants or contracts, 
and to test the durability of aircraft gas turbine engine technology in the 
SSME environment. 


EXISTING TECHNOLOGY 

Through a series of NASA contracts with Pratt & Whitney Aircraft, the 
thin-film temperature senscr technology was developed for aircraft gas turbine 
engines (refs. 2 to 4). In this technology the turbine blades coated with a 
proprietary NiCoCrAlY anticorrosion coating are oxidized to grow an adherent, 
electrically insulating surface film of aluminum oxide. This film is augmented 
with sputtered AI 2 O 3 and then the thermocouple is sputter deposited (fig. 2 ). 
Junctions between the thin films and the lead wires are formed by hot- 
compression diffusion bonding. These thin-film sensors have performed success^ 
fully in engine tests to 1800 for over 60 hr with less than 0.02 percent 
drift per hour. 

STATUS OF WORK 

An initial goal is to augment the in-house capability in thin-film sensor 
technology. A thin-film sensor laboratory is being installed in a refurbished 
clean room, and new sputtering and photoresist exposure equipment is being 
acquired. Also, the procurement of material samples is in progress, and an 
effort is under way to put the initial coating on the samples. The MAR-M 200 
(-i-Hf) and MAR-M 246 (+Hf) along with other materials are to be coated with 
NiCoCrAlY and NiCrAlY coatings. 

In an on-going research activity through a grant with UCLA, various coat- 
ings and their insulating films are being investigated for use in sensor devel- 
opme*^t. An impetus of this activity is to develop a generic coating material 
tha s comparable to the proprietary MCrAlY coating. Some coating deposition 
methods such as activated reactive evaporation and chemical vapor deposition 
are being explored. The emphasis is on improving the processing techniques for 
greater processing yield. 
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Another activity is the testing of existing thin-film thermocouple tech- 
nology in an SSME environment. Test blades will be made from MAR-M 200 (+Hf), 
for which the gas turbine type of thin-film sensor was developed. These blades 
are to be instrumented with two thin-film thermocouples deposited on one side 
of the airfoil and a wire thermocouple embedded from the other side of the 
blade (fig. 3). These instrumented blades will be tested in the turbine blade 
tester at the Marshall Space Flight Center (fig. 4). Besides testing the 
survivability of the aircraft turoine type of senscr in an SSME environment, 
this program v/ill provide a measurement of the cyclic thermal gradient imposed 
on the blades. 


FUTURE PLANS 

Presently the SSME turbine blades are coated only with NiCrAlY bond coat- 
ing, which provides minimal thermal protection. As more information on struc- 
tural integrity and durability becomes available, the blade and the coating 
material may change. In anticipation of future improveineats many varieties of 
alloys and coatings will be tested for compatibility with thln-film sensors. 

The thin-film technology extends to other sensors such as the strain 
gauge. Figure 5 shows a thin-film strain gauge sputtered onto an aluminum 
oxide beam. This thin-film strain gauge will be tested in a high-temperature 
oven to determine the performance characteristics and adhesion properties of 
the sensor deposited on a ceramic surface. 

REFERENCES 

1. Holmes, R.R: Vacuum Plasma Coatings for Turbine Blades. Proceedings of 
MSFC Advanced High Pressure O 2 /H 2 Technology Conference. June 1984. 

2. Grant, H.P.; and Przybyszewski , J.S.: Thin Film Temperature Sensor. NASA 

CR-159782. 1980. 

3. Grant, H.P.; Przybyszewski, J.S.; and Claing, R.G.: Turbine Blade Tempera- 
ture Measurements Using Thin Film Temperature Sensors. NASA CR-165201, 
1981. 

4. Grant, H.P.; Przybyszewski, J.S.; Claing, R.G.; and Anderson, W.L.; Thin 
Film Temperature Sensors, Phase III. NASA CR-165476, 1982. 
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THE STABLE ADHERENT AI 2 O 3 INSULATING UYER IS OBTAINED 
BY AT LEAST 50-hr OXIDATION (AT 1300 K) OF ^HE COATING. 
FOLLOWED BY AI 2 O 3 SPUHERING. 


Figure 2.- Thin-film temperature sensor on turbine blade. 



FLANGE 


(a) Blade holder assembly. Th>'ee blades, 
each blade tc have two thin-him thermo- 
couples and one embedded-wire thermo- 
couple. 



(b) Potential location. Twothin- 
hlm thermocouples on pres* 
sure or suction side, with wire 
thermocouple eiitbedded from 
opposite side. 

Figure 3,- Conceptual sketch of blade holder assembly. 
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W85-27969 

OPTICAL STKAIH HEASURBMBHT SYSTEM DBVBLOPMBBT* 


Christian T. Lant 
Sverdrup Technology, Inc 
Middleburg Heights, Ohio 44130 


Investigations of physical phenomena affecting the darability of S3MB 
con^onents require measurement systems operational in hostile environments. 

The need for such instrumentation bore this contract to define and put into 
operation an optical strain measurement system. Sverdrtip Techno] ogy, Inc., in 
cooperation with the Lewis instrument research laboratory, is currently devel- 
oping this optical strain measurement system based on the speckle-shift method 
of Ichirou Yamaguchi fThe Institute of Physical and Chemical Research, Japan). 
This is a noncontact, automatic method of measuring surface strain in one 
dimension that corrects for error due to rigid-body i^otion. It provides a 
gauge length of 1 to 2 mm and allows the region of inte**est on the test speci- 
men to be mapped point by point. The output will be a graphics map of the 
points inspected on the specis^n; data points will be stored in quasi-real 
time. This task is the first phase of a multiphase effort in optical strain 
measurement. Phase 2 will extend the technique to measure strain components 
in two dimensions. 


BASIC PRINCIPLES 

Theory shows that the difference between the speckle displacements for a 
pair of synanetrically incident laser beams is directly proportional to the 
surface strain parallel to the plane of incidence and independent of the 
translational and rotational components of motion. The speckle pattern created 
by the test specimen is interpreted as high-order interference fringes result- 
ing from a random diffraction grating, being the natural surface roughness of 
the specimen. Strain r.nduced on the specimen causes a change in spacing of 
the surface roughness, which in turn shifts the position of the interferOi.ce 
pattern (speckles). 


EXPERIMENTAL CONFIGURATION 

In the experimental setup an argon-ion laser and a beam splitter provide 
the two symmetrical beams, which will each independently expose a linear photo- 
diode array before and after stress is applied to the specimen. The ''before*' 
and "after" exposure patterns of the speckle in the array are cross-correlated 
fot each beam to give the spatial shifts. After the shift from each beam is 
computed, the difference is taken to cancel out components due to rigid-body 
motion; the remaining shift is due to surface strain. The sensitivity of this 
strain gauge depends on the ratio of the pitch of the diodes in the linear 
array to the distance from the array to the test specimen. 

The test specimen will be mounted on a stress fatigue testing machine, 
with specimen dimensions of approximately 100 mm x 20 mm x 1 mm. The center- 
to-center spacing of the diodes will be 15 |im with a sensor- to-specimen 
spacing of the order of a meter. Thus the resolution of the setup will be 10 


*Work performed under Lewis task number 83-24-01. 
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to 20 microstrains, with an unlimited range because the increatental speckle 
displacements can be summed by using successive exposures. The system %fill be 
used for mechanically and thermally induced strain, and the experimentation 
will cover a temperature range of 20 to 600 deg Celsius. Testing will deter- 
mine the amount of decorrelation of the speckle pattern at high tessera tures . 
The upper tei^erature limit will be a function of this decorrelation. A gauge 
length of 1 to 2 USB corresponds to the spot size of the laser beaais. The 
region of interest can be mapped by relocating the laser spots. Utilization 
of the natural surface structure of the specimen eliminates any need to arti- 
ficially prepare the surface of the specimen before testing. However, for 
high-tenq>erature strain measnirements , cycling through the teiQ>erature range 
may be necessary to stabilize oxidation of the surface. 

OPERATIOII 

The system will be designed for a minimal amount of operator input; the 
procedure will be run by a microcomputer system. The operator will be required 
to position the laser spots through the terminal of the coaqmter and to siain- 
tain this position with the aid of an area photodiode array. This area array 
will provide a graphics display of the test specimen and laser spots for a CRT 
and store location of the data points taken on the specimen. Motorized mirror 
mounts will be used to direct the beams onto the test specimen. This combina- 
tion will ensure that relative movement between the test apparatus and the 
laser does not influence the accuracy of the strain measurement. 

REFERENCE 

1. Yansaguchi, I.: A Laser-Speckle Strain Gauge. J. Phys. E: Sci. Instrum. , 

vol. 14, 1981. Printed in Great Britain. 
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SCOPE 


Multi-task effort: 


This task: One dimensional strain map 


/uture task: Measure strain of SSME components 

In two dimensions 


TASK OBJECTIVES 


Operationalize Yamaguchis* technique 


- one dimensional strain measurement 


- real or quasi-real time output of strain 
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Correlation Technique 




maxim* ti of Rs^valaa of ax 

fCxIegCxvaxI 
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STRAIN IS COMPUTED 
FROM CORRELATION RESULTS : 

1. Cross-correlation gives speckle displacement AX 
d. Strain is directly proportional to AX : 

e= -.AX_ 

2L8in9 
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HSAT FLUX SENSOR CALIBRATOR 


N85-27970 


Curt H. Liebert 

National Aeronautics and Space Administration 
Lewis Research C«iter 
Cleveland, Ohio 44135 


The heat flux to space shuttle main engine (SSMB) turbopuiq> turbine blades 
may be as high as 10^ U/m^. This heat flux is 50 to 100 times that encoun- 
tered in aircraft engines. The heat flux causes thermal transients that are of 
the order of 1 sec as temperature varies from perhaps 1500 K to 100 K. It is 
suspected that these transients cause durability problems in the turbine 
blades. To quantitatively evaluate the effect of these transients, heat flux 
sensors or gauges are needed to obtain data to verify analytical models 
(ref. 1). Because of these very high heat fluxes a new system must be devel- 
oped for the calibration of heat flux sensors mounted on these blades. 

OBJECTIVE 

The objective of this effort is to design and fabricate a system for 
steady-state and transient calibration and durability testing of heat flux 
sensors for use in SSHE turbine blades. Development will rely on experiment 
and analysis from many interdisciplinary fields of research including thln-film 
and wire thermocouples, materials, heat transfer, optics, and electric arc 
phenomena. Building on the techniques available in these fields, calibration 
and durability testing skills will be developed for measuring temperature and 
heat flux on SSMB turbines. 


LITERATURE SEARCH 

The literature presents siany methods for heat flux calibration. Most of 
the methods are based on the use of thermal radiation from solids; other meth- 
ods are based on conduction or convection heat transfer modes (ref. 2). All 
of these methods are useful for calibration only at heat flux levels of one or 
two orders of magnitude less than those estimated for SSME turbopump blades. 
However, new developments in industrial electric arc engineering in the past 
several years have made electric arcs useful as a radiative heat source for 
calibration at the higher heat fluxes generated in SSME turbines. This devel- 
opment (fig. 1) is based on a patented electric arc lamp design (fig. 2 , 
ref. 3) that produces optical radiation focused with a truncated elliptical 
mirror (ref. 4). The electric arc lamp design was based in p<.j‘t on research 
done at the Lewis Research Center during the 1960's and 1970* s (ref. 5). 
Adjustable irradiances can be focused on the entire area of a pressure or suc- 
tion surface of an SSMB turbine blade at levels to 10^ W/m^ with a uniformity 
of 1 to 2 percent. The lamp will be remotely operated in either a pulsed or 
steady-state mode. Its performance is advertised as repeatable and indepen- 
dent of the operator or of the sensor being calibrated. Lamp life is several 
hundred hours. 


CALIBRATOR DESIGN AND INSTALLATION 

The calibrator will be used in a calibration and durability testing facil- 
ity that includes equipment shown in figure 3. The calibrator consists of 
(1) the arc lamp, (2) a high-speed positioning table for placing standard and 
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special heat flux sensors in the incident beam of radiant heat flux, (3) pro- 
vision for cooling the blade and special sensors inserted in the blade, (4) a 
computer for controlling the positioning table (local and remote), storing 
electrical output values from sensors, and calculating heat flux from these 
values, and (5) a pyrometer for measuring sensor surface teng>eratures. 

CALIBSATIOH OF THE CALIBRATOR 

The intensity of the arc will be calibrated with commercially available 
heat flux sensors that have been factory calibrated by techniques traceable to 
those used by the Rational Bureau of Standards. 

DBHONSTRATIOH TESTS 

Facility tests will include measurement of sensor electrical output as a 
function of known radiant heat flux input. Both transient and steady-state 
tests will be performed. 

Heat transfer analysis and thermal stress experiments are proceeding for 
the design and fabrication of, and the reduction of data from, special heat 
flux sensors mounted in turbopunq> airfoils. The special sensors will be used 
with the calibrator and other equipment to demonstrate the operational aspects 
of the facility. 
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FEASIBILITY OF MAPPING VELOCITY FLOW FIELDS IN SSME 

POWEEHEAO BY LASER ANENOMETRY TECHNIQUES* 

D.G. Pclacclo, L.K. Sharma, T.V. Ferguson, 

J.M. Haram, and S. Plnkowski 
Rockwell International 
Canoga Park, California 91304 


The requirements of higher performance for the space shuttle main engine 
(SSME) and future Earth-to-orbit propulsion systems dictate the use of effi- 
cient liquid- fueled rocket engines with cycles that require lightweight, high- 
pressure, conq>act powerhead configurations. These couf ig**rations typically 
exhibit conq>lex, nonuniform hot-gas internal fluws. Su':h propulsion systems 
must be designed with an emphasis on long life, higt’ lability, and minimum 
maintenance. In-depth characterization and inderrt. , of the powe..head flow 
field is essential for optimal design of this key ooraponent . bASA is 

currently planning hot-fire engine tests for the S£H, Enhancement Tech- 

nology Program to establish a detailed data base from which to address power- 
head flow issues and to verify advanced, three-dimensional, fluid-dynamic 
analysis models currently under development. 

Because of the flow environment associated with the SSME powerhead pres- 
sure (3000 psia), temperature (IPOO **R) , and mechanical complexity and the 
high-vibration test stand environment, detailed flow measurements are at best 
difficult to make. During the past year the NASA Lewis Research Center initi- 
ated a study with Rocketydne that addresses the feasibility of using laser 
anemometry techniques to map velocity flow fields in an SSME powerhead. This 
is in support of the proposed model verification experiments. 

In the on-going study three engine powerhead conq>onent flow environments - 
(1) the high-pressure fuel turbopump prebumer, (2) the fuel turbopump turbine 
rotor and stator region, and (3) the 180* turnaround duct - are being consid- 
ered. Flow parameters to be measured by the anemometry techniques are time- 
averaged values of the velocity magnitude and flow direction, turbulence 
intensity, velocity component correlations, integral time scale, and tart fierce 
spectrum. 

Key technical issues such as identifying feasible means of optical access 
to the high-pressure, high-temperature measurement flow regions, and measure- 
ment system conq>atibi1 ity with the test environment are addressed. At the 
conclusion of the study a measurement system is to be defined that will include 
layout sketches and specifications of component hardware that will be compat- 
ible with mapping the flow field in one r.^ the regions of interest. 


*Work performed under NASA contract NAS3-24356. 
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